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ABSTRACT 

Context. The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. Studying 
the resolved chemical substructures of this complex region provides important insight into the chemistry of high-mass star-forming 
regions (HMSFRs), as it relates to their evolutionary states. 

Aims. The goal of this work is to resolve the molecular line emission from individual substructures of Orion-KL at high spectral and 
spatial resolution and to infer the chemical properties of the associated gas. 

Methods. We present a line survey of Orion-KL obtained from combined Submillimeter Array (SMA) interferometric and IRAM 
30 m single-dish observations. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species 
(25 isotopologues), including 11 complex organic molecules (COMs). Spectra are extracted from individual substructures and the 
intensity-integrated distribution map for each species is provided. We then estimate the rotation temperature for each substructure, 
along with their molecular column densities and abundances. 

Results. For the first time, we complement 1.3 mm interferometric data with single-dish observations of the Orion-KL region and 
study small-scale chemical variations in this region. (1) We resolve continuum substructures on ~3"angular scale. (2) We identify 
lines from the low-abundance COMs CH 3 COCH 3 and CH 3 CH 2 OH, as well as tentatively detect CH 3 CHO and long carbon-chain 
molecules C^H and HC 7 N. (3) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., 
nitrogen-bearing species) or the compact ridge (e.g., oxygen-bearing species like HCOOCH 3 and CH 3 OCH 3 ), the distributions of 
others do not follow this segregated structure (e.g., CH 3 CH 2 OH, CH 3 OH, CH 3 COCH 3 ). (4) We find a second velocity component of 
HNCO, SO 2 , ^^S 02 , and SO lines, which may be associated with a strong shock event in the low-velocity outfiow. (5) Temperatures 
and molecular abundances show large gradients between central condensations and the outfiow regions, illustrating a transition be¬ 
tween hot molecular core and shock-chemistry dominated regimes. 

Conclusions. Our observations of spatially resolved abundance variations in Orion-KL provide the nearest reference source for hot 
molecular core and outfiow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs. 

Key words. Stars: formation; Stars: massive; ISM: lines and bands; ISM: molecules; ISM: abundance; Submillimeter: ISM 


1. Introduction 


The birth sites of high-mass star-forming regions (HMSFRs) 
are highly complex structures, consisting of different gas and 
dust cores, embedded in a less dense gas envelope. While 
the complexity of dynamic processes within these regions 
has been characterized well (e.g., [Zinnecker & Yorke||2007[ 
[Beuther et al.||2007bt |Tan et al" 2QI4| ), our knowledge of the 
chemical properties of these regions is still poor. Since specific 
molecules can act as tracers of phy s ical processes (e.g.,|Schilke 
et al.l TW7}JChamley et al.||I997[ [Bachiller et al. |20()It |Viti 


et al.|2004[|Garrod||20I3| ), characterizing the chemistry of gas 


associated with the early stages of high-mass stars is crucial for 
understanding how these stars might form. 


Line surveys of HMSFRs in mm/sub-mm at mospheric 
windows exhibit rich molecular emission spectra (see |Herbst &| 
van Dishoeck|2009] for a comprehensive review). The multitude 
of detected lines allow for an unbiased look at the chemical 
composition of a particular region, from which its physical 
properties can be inferred, including temperature and density. 
However, HMSFRs are typically distant (several kpc), so that 
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spatially disentangling the physical and chemical substructures 
is challenging with the resolution of most facilities. This makes 
the Becklin-Neugebauer and Kleinmann-Low Nebula in Orion 


(Orion BN/KL), at a dist ance of 414 + 7 pc ( [Sandstrom et al. 
2007[|Menten et al.|2007| ), ideal for studies of the physical and 


chemical properties of HMSFRs. As the nearest HMSFR, it has 
been extensively studied at both far-infra red to (sub) millimetre 
wavelengths accessible from space (e.g., Olofsson et al.||2007} 


|Bergin et al. |20 1 0[[Crockett et al.|20I4| ) and at radio, millimetre, 
and submillimetre wavelengths accessible from the ground. 
A broad range of studies on this source have been conducted 


Sutton et al.||I985t Blake et al. |I987t |Turner||I989l GreavesI 

& White||I99It ISchilke et al.||I9971 |200It |Lee & Cho 

12002 

Comito et al. 2005[ |Remijan et al. 2008[ Tercero et al. 

(2010) 

or interferometers (e.g., Wright et al. 1 19961 |Blake et al. 

11996 

Beuther et al.|2005| 2006[ |Friedel & Snyder|2008l Guelin et al. 

2008|). Strong emission lines in these observations have revealec 

complex dynamics (e.g., |Sutton et al.||I985 

|Blake et al. 

1987 

|Turner|I989t|Beuther & Nissen|2008l|Beut 

her et al.|20I0|) and 

chemistry (e.g., Wright et al.||I996l Blake et al.||I996 E 

leuther 

let al.||2005 20061 Friedel & Snyder 20081 

Guelin et al.||20081 

[Widicus Weaver & Friedel||20I2l |Friedel & Widicus Weaver| 
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I2OT2I ). 


Both single-dish-only and interferometer-only observations, 
however, have limitations. Single-dish observations have low 
spatial resolution, making it difficult to detect spatial variations 
of chemistry on small scales. In contrast, interferometric 
observations can probe spatial variations on small scales, but 
suffer from filtering-out to large scale structures (the “missing 
short spacing” problem; Bajaja & van Albada 1979 ), making 
measurements of the chemistry of more extended features unre¬ 
liable. In this paper, we aim to study the chemistry of Orion-KL 
by obtaining a combined dataset from an interferometer 
(Submillimetre Array, SMA) and a single-dish (IRAM 30 m) 
telescope. These combined data, which resolve small-scale 
substructures and are sensitive to more extended components, 
allow us to accurately map the distribution of molecular species, 
determine their abundances, and ultimately analyse spatial 
variations in the chemistry. 


One unique feature of HMSFRs is hot molecular cores, 
which are compact objects (< 0.05 pc) that e xhibit high temper - 
atures (> 100 K) and densities (> 10^ cm~^; ^au et al.[|2011 1 . 


Numerous spectral surveys (see, e.g., |Turner| |1989[ |Schilke 
et al.|200n|Belloche et al.|20'09l|Herbst & van Dishoeck|20091 


Tercero et al.||2010| ) have revealed that hot molecular cores are 

characterized by a large abundance of simple organic molecules, 
including H 2 CO (formaldehyde), HCOOH (formic acid), and 
complex organic molecules (COMs, organics containing > 6 
atoms), such as CH 3 OH (methanol), HCOOCH 3 (methyl 
formate), and CH 3 OCH 3 (dimethyl ether). Although some 
larger COMs have lo w abundances in ho t molecular cores, e.g., 
CH 3 CH 2 OH (ethan ol; [Rizzo et al.|2001| ), CH 3 COCH 3 (acetone; 
Snyder et al.|20Q2| ), CH 3 CHO (aceta ldehyde; [Nummelin et ak] 


2000D, (CH 2 bH )2 (ethylene glycol; [Hollis et al.||2002|), and 
CH 2 OHCHO (glycolaldehyde; [Hollis et al.[[2000r|2004f , they 
have been previously detected towards the Galactic-centre 
hot molecular core source SgrB2 (N-LMH). Because of their 
high degree of chemical and structural complexity, the precise 
origins and formation mechanisms of these COMs is still a 
subject of debate. Nevertheless, understanding these topics is 
particularly important for understanding the chemical evolution 


of the interstellar medium (ISM; e.g., Herbst & van Dishoeck 
[2009[). Recently, a number of interferometric observations 


of COMs in Orion-KL (e.g., Eriedel et al. 

20051 [Eriedel & 

Snyder[[2008[ Guelin et al.|2008[[Pavre et al. 

2011a|b| Widicus 

Weaver & Eriedel 2012) have shown spatia 

-related chemical 


differentiation. However, the relatively low sensitivity of these 
observations means that the lowest-abundance COMs (e.g., 
CH 2 OHCHO, CH 3 CH 2 OH, and CH 3 CHO) are below the 3o- 
detection limit or are blended by strong neighbouring lines. 


In addition to the hot molecular core chemistry, the shock 
chemistry driven by outflows is another important aspect of 
chemistry in HMSFRs. The region in Orion-KL that is most 
affected by outflow and shock events is called the “plateau” 
( Blake et al.[[19p[ [Wright et al.[p^^ [Lerate et al.[[2008| 
|Esplugues et al.[[2013[ ) and is characterized by spectral lines 
with broad line wings. Two roughly perpendicular outflows are 
seen in Orion-KL, traced by the proper motion of H 2 O masers 
at 22 GHz ( [Genzel et al.|1981[ ). One is a high-velocity outflow 
(30-100 kms“\ along the north-west-south-east (NW-SE) 
axis); this outflow is also detected in OH (hydroxyl) masers 
( Norris[[ 19841 [Cohen et al.[[2()06|), broad bipolar CO line wings 


( Zuckerman et al.|1976t[Kwan & Scoville|1976|[Erickson et al. 


1982, Masson et al.[[1987 

Chemin & Wright[[1996[ [Beuther & 

Nissen|2008t 

Zapata et al. 

2009), and lobes of poorly collimated 

shock-excitec 

H 2 ([Beckwith et al.[[1978[ [Scoville et al.[[1982l 

Taylor et al. 

19841 Allen & Burton 

19931 Sugai et al. 

[1994). 


Although the source responsible for driving the outflow is still 
uncertain, it has been suggested that the outflow is only about 


500-700 years old ([Lee & Burton|2000t[Doi et al.|2002 

: Gomez 

et al. 20051 Nissen et al.[[2007t [Gomez et al. 

2008 

Beuther 

& Nissen 120081 Bally et al.[[2011l [Goddi et al. 

bout 

) [Nissen 

et al.[[2012[). The other is a low-velocity outflow (~ ] 

8 kms“\ 


aV ~ 35 kms ^ along the north-east-south-west (NE-SW) 


axis; [Genzel et al.[[1981t [Gaume et al.[p^98] ), whi ch is also 
detected in SO (Plambeck et al.[ 1982|), SiO masers ([Wright & 
[Plambeck][1983t Plambeck et al.[ 20 ^ [Greenhill et al.[[ 20 T^ 
and H 2 O maser ( [Greenhill et al.|2013[ ). This outflow is believed 
to be closely associated with an expanding gas shell or torus 
around Source I, but whether its age is > 3000 years ([Menten 


& Reid| 1995 [[Greenhill et al.|1998t[Doeleman et al.| 1999 [[Reid 


et al.|2607 ) or much younger ( Plambeck et al.[[200^ Greenhill 

et al.|2013 ) is still not clear. Another goal of this study is to look 
for additional chemical signposts of these outflows, apart from 
the known masers, in order to better understand their physical 
and chemical properties. 


In this paper, we present the first combined single-dish and 
interferometric line survey at 1.3 mm towards Orion-KL and 
study the chemical properties of this region at a linear resolution 
of ~1 200 AU (at a distance of -414 pc). In Section we de¬ 
scribe the observations and their combination. In Section we 
present a continuum map and maps of the distribution of differ¬ 
ent molecules, including simple molecules and COMs, where the 
presences of COMs are identified by synthetic spectrum fittings. 
The spatial distributions, column densities, and abundances for 
all the detected species are derived in Section]^ A discussion of 
the chemistry in the hot molecular core and outflows of Orion- 
KL can be found in Section [5l 


2. Observations and data reduction 


2.1. Submillimetre Array (SMA) 

We carried out observations of Orion-KL with the SMA on Eeb. 
15, 2005 and Eeb. 19, 2005. At 1.3 mm (220/230 GHz) in the 
compact configuration with seven antennas, the baselines ranged 
between 16 and 69 m. The short baseline cutoff causes source 
structures > 20"to be filtered out. The primary beam is 52", and 
the phase centre of the observations is 0^2000 = 05^35”^ 14.5^ and 
^2000=-5°22'30.45" with Visr ~ 7 kms“^. Bandpass calibration 
was done with Callisto, and the flux calibration is estimated 
to be accurate within 15%, based on the SMA monitoring of 
quasars. Phase and amplitude calibration was performed via 
frequent observations of the quasar 0607-157. The zenith opaci¬ 
ties, measured with the NRAO tipping radiometre located at the 
Caltech Submillimetre Observatory, were excellent during both 
tracks with t(L3 mm) - 0.03 - 0.04, allo wing us to observe a t 
both 1.3 mm and 440 jum simultaneously ( [Beuther et al.|2006] ). 
More technical descrip tions of the SM A and its calibration 
schemes can be found in Ho et al. ( [2004| ). The receiver operated 
in a double-side band mode with an intermediate frequency 
of 4-6 GHz so that the upper and lower side bands were 
separated by 10 GHz. The correlator had a bandwidth of 2 GHz, 
and the channel separation was 0.8125 MHz, corresponding 
to a velocity resolution of ~ 1.2 kms“^. Double-sided band 
system temperatures corrected for the atmosphere were be- 
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tween 150 and 300 K (depending on the elevation of the source). 


The initial flagging and calibration was done with the IDL 
superset MIR, which is originally developed for the Owens 
Valley Radio Observatory and adapted for the SM^Apl The 


imaging and data analysis was conducted in MIRIAD ^Sault| 


et al.||l995| ) and CASA0 Between the observations 


on two 

nights, the frequency setups were shifted by 0.4 GHz, yield¬ 
ing a slightly larger total bandwidth for the combined dataset 
(218.870-221.234 GHz for the lower sideband, and 228.870- 
231.234 GHz for the upper sideband). The synthesised beams of 
the combined data are 3.49" x 2.57" for the lower sideband and 
3.68" X 2.69" for the upper sideband. 


2.2. Single-dish observations with the IRAM 30 nn telescope 

In addition to the high spatial resolution data from the SMA, 
we mapped an area of 1' x 1', with the IRAM 30 m telescope, 
to complement the short-spacing information missing from 
the interferometric observations. Observations were performed 
in the on-the-fly mode on Nov. 30, 2012, covering a broad 
bandpass (8 GHz bandwidth for each sideband) of EMIR 
at 220/230 GHz. Observations were conducted under good 
weather conditions (Tsys = 266 K, r ~ 0.152 at 220 GHz). The 
focus was on Uranus, and pointing was checked on Uranus 
and 0420-014. Using a forward efficiency of 92% and a main 
beam efficiency of 58%, we converted the data from antenna 
temperature (K) to flux (Jy) by using an efficienc}0of 7.8. The 
beam of the 30 m telescope at 1.3 mm is ~ 11.3" in the upper 
sideband and ~ 11.8" in the lower sideband. The observations 
cover frequencies 213.370-221.150 GHz and 229.051-236.831, 
with a velocity resolution of 0.255 kms“^. We used the Gilda^ 
software for data reduction, and then converted the 30 m data to 
the MIRIAD data format. 


2.3. Data merging 

The 30 m single-dish spectral-line data were first converted 
to visibilities, and then combined with the SMA data us¬ 
ing the MIRIAD package task UVMODEl]^ To recover 
large-scale emission, while simultaneously maintaing high 
spatial resolutions for resolving small scale structure, we 
chose an intermediate weighting between natural and uniform 
(“robust”= -0.5 value in MIRIAD). However, as continuum 
information is not available from the 30 m observations, only 
the spectral-line data can be corrected in this way. 

The combined data have a field of view of 52" and a beam 
of 5.55" x4.35" with position angle P.A. = -i-8° in the lower 
sideband and 5.44" x 4.43" with P.A. = -\-2° in the upper 
sideband. The frequency overlap between the SMA and 30 m 
data results in the combined data having a frequency range 
218.870-221.150GHz and 229.051-231.234GHz. Eigure Q 
presents an amplitude-versus-projected uv distance plot, show¬ 
ing the good combination of the single-dish and interferometer 
data. Missing projected baselines shortwards of 10 kT from the 

^ http: //cfa-WWW. harvard. edu/~cqi/mircook.html 
^ http://casa.nrao.edu 

^ http://www.iram.es/IRAMES/mainWiki/Tram30mEfficiencies 
4 http://www.iram.fr/IRAMFR/GILDAS 

^ https://www.cfa.harvard.edu/sma/miriad/manuals/SMAuguide/ 
smauserhtml/uvmodel. html 




Fig. 1: Amplitude in Jy versus the projected baseline showing 
the combination of the interferometer and single-dish data, from 
the upper sideband (top) to lower sideband (bottom). Red and 
blue dots present data from the SMA on Feb. 15, 2005 and Feb. 
19, 2005, respectively, showing the lack of baselines shorter than 
lOk/1. Black dots are data converted from 30 m single-dish ob¬ 
servations, which fill in the missing short spacing information. 


SMA have been complemented by using the 30 m data. After 
combination, the Icr rms per 1.2 kms“^ channel is measured to 
be 80 mJy beam“^ for the lower sideband and 70 mJy beam“^ 
for the upper sideband. 


3. Observational results 

3.1. Continuum emission 


Orion-KL has a rich spectrum of molecular emission fines, so 
extracting the continuum requires some care. To construct our 
continuum map, we averaged the line-free parts of the upper 
and lower sidebands from the SMA-only data. In Figure we 
present the resulting 1.3 mm continuum image as contou rs over¬ 
laid on a m ap of the outflow traced by shock-excited H 2 ( [Nissen 
et al.|2QQ7] ). We detected the previously known compact contin¬ 
uum substructures of this region and moreover, resolved some of 
these substructures into multiple distinct condensations. We fist 
the following resolved condensations in Table 


- The brightest hot core (HC): Located close to a well-studied 
infrared source IRc2 ( [Downes et al. 1981[ [Wynn- Williams 


et al.|198^ , the brightest HC is the first clump in Orion-KL 
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detected by the high spatial resolution mm observations |et al.|20lTa| ), it has not been characterized in detail. 


([Masson et al.|1985}|Mundy et al.|1986{ [Woody et al.|198^ 

[Wright et al.|1992[ )7 and it exhibits the highest peak specific 

intensity (Jy beam“^) in the 1.3 mm SMA continuum map. 
With the temperature estimated to be 100-300 K (e.g., 
[Wilner et al.|T994l [Wilson et al.|2000[[Beuther et al.|2005l ), 

this source has been found to be abundant in hydrogenated 
species and nitrogen- (N-) bearing molecules ([Sutton et al. 


1995 


Blake et ^[1996[ [Wright et al.[[1996t [Wilson et al. 
2000[). With even higher resolution (0.7", ~300 AU), obser¬ 


vations from SMA at 865 yum/348 GHz expose additional 
substructures within the HC ( [Beuther et al.|2004] ). Presented 
as black contours overlaid on the 1.3 mm colourmap in the 
inner panel of Figure]^ these 865 //m structures are: SMAl 
(a high-mass protostellar source, sugges ted to be the driving 
source of the high-velocity outflow; Beuther & Nissen 
[2008[ ), Source I (located at the centre of the SiO masers 
and believed to drive the bipolar low-velo city outfiow; 
[Menten & Reid[p^95l [Matthews et al.[[MT 0 | ), Source N (a 
Herbig Ae/Be or mid-B star with a luminosity of ~ 200QLq; 
Menten & Reid[[1995i [Greenhill et al.][2004[ [Gomez et al. 


2005 


et al. 


[Rodriguez et al.[[2005t [Goddi et al.[[2011bt [Nissen 
2012[ ), and hotcore. Hereafter, we use “HC” to denote 


the whole region resolved by 1.3 mm continuum at a spatial 
resolution of ~1 200 AU, the “hotcore” to denote the posi¬ 
tion resolved by 865 yum continuum at a spatial resolution 
of -300 AU, and the “hot molecular core” to describe an 
evolutionary status. 

mm2: Located - 7" west of the HC, and - 5" south of the 
BN object. Coincident with infrared sources IRc3, IRc 6 , 
and IRc20 ( [Dougados et al.|19'93 ), t his condensation is also 
known as the western clump (WC; Wright et al.[[1992 ), or 
the north-west clump (NW C; [Blake et al.[[ 19961 [Tang et al "' 


2010 Favre et aTl[2011a ). The mass (2-3 Mq), H 2 column 


density, and high temperature (a few hundred K) in this 
source suggest that it may be the site of ongoing low-mass 
star formation ( [Schreyer & Henning|1999[ ). 

mm3 a and mm3b: These two condensations are located 
south-west of the HC . Coincident with infrared IRc5 
( Dougados et al.[[l993[ ), they were previously considered 
to be a single clump called the south-west clump (SWC; 
[Tang et ^ 2010[ ), or Cnt C ([Murata et aL][1992[), situated 
at the compact ridge (C R; [Blake et al.[[1996[ [Wright et al. 


[ 19961 [Thng et al.|2010[ ). Similar dual continuum peaks have 
been detected in other high spatial resolution observations 
of HMSFRs (jTang et al.|2010t [Favre et al.|2011at[Friedel &[ 


Widicus Weaver|2Qll[ ), so we^consider them to be separated. 

The CR is known to have different kinematics and ch emistry 


The north-east clump (NE): a strong continuum emission 
peak which is adjacent the HC and the north extended ridge 
(ER; [Sutton et~n^[1985[ [Blake et ar|[1987[ ). This clump is 
usually filtered out in high spatial resolution observations, 
but can be clearly detected in interfer ometric observations 
made with large scale configurations ( [Murata et al. [1992} 
Tang et al.|2010[ ). Comparison with the existing single-dish 


maps suggests it is part of a large-scale dusty filament along 
the NE-SW direction (e.g., detected by JCMT at 450 yum, 
Wright et al.|[T992l SC UBA at 450 jum and 850 ixm, |m 


Francesco et al.||2608|l. In addition, this clump is in the 


plateau of the low-velocity outfiow, so its chemistry may be 
affected by shocks (Section [53] ). 


Because of the limited UV coverage and field of view, the 
high-velocity outflow and the BN object are not detected above 
5(T in the SMA-only continuum map. To study the chemistry in 
these regions from the combined data, we take the positions of 
the high-velocity platea u (Visr ~ 11 kms“^) and BN object from 
[Beckwith et al. ( [1978[ ) and denote them as OEIN, OEIS, and 
BN in Eigure[^ 


than in the HC: it is a vast regi on (~ 15"-25"; Plume[ 
[et al.[[20f^ [Esplugues et aL|[2013[ ) with lower temperature 
( [Genzel & Stutzkif 1989[), is dominated by oxygen- (Q-) 


bearing molecu les ( Blake et"ar][1987[ Wright et al.[[1996( 
[Liu et al.|2002[ ), and is believed to be the region where two 
outflows in t eract with the d ense ambient material ( [Blake[ 
[et al.[[198T| [Liu et al.][2002[ ). Therefore, to ensure that we 
fully characterise its chemical properties, the following 
study focuses separately on mm3 a and mm3b, instead of 
considering CR as a whole. 

The southern ridge (SR): the southernmost condensation 
detected in the 1.3 mm continuu m map. Although it has been 
detected in other observations ( [Eisner et al.[[2008] [Eavre[ 


Table lists the nominal absolute positions of these sub¬ 
structures, as well as the peak specific intensity (Jy beam“^) 
for each continuum substructure, obtained from the SMA-only 
1.3 mm continuu m map (Eigure for more detail about this 
table, see Section 4.2.2[ ). The HC is the brightest condensation; 
mm2, mm3a, and mm3b are fainter; while NE, SR, OEIN, and 
OEIS are even weaker. 


3.2. Spectral line emission 
3.2.1. Line identification 

After imaging the whole SMA-30m combined data cube, we 
extracted the double sideband spectrum from each substructure 
denoted in Eigure and show them in Eigure The majority 
of lines are strongest towards the HC, so we use the spectrum 
from the HC to identify lines. However, line identification may 
be ambiguous in a region such as the Orion-KL HC, because 
of a high density of lines (that can have different line widths 
according the species), also because the - 1.2 kms“^ spectral 
resolution of the observations does not allow us to separate 
all the blending. Eor strong lines, we first compare their rest 
frequencies with an existing single-dish survey catalogue of 
Orion-KL by [Sutton et al.[ ( [1985[ ). Next, we cross-check their 
rest frequencies with the molecular database at “Splatalogue’r 

in i'. , ^ 


(a compilation of the Jet Propulsion Laboratory (JPL ' Pickett 


[et al.[p^98] ), Cologn e Database for Mole cular Spectroscopy 
catalogues (CDM ^ ’ 


Muller et al. 


2005), and Lovas/NIST 


catalogues, [Lovas 2004), and confirm the identification if the 


o www.splatalogue.net 
^ http://spec.jpl.nasa.gov 
^ http://www.astro.uni-koeln.de/cdms/catalog 
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Fig. 2: Substructures resolved by SMA-only observations at a spatial resolution of ~ 1200AU. White contours show the 1.3 mm 
continuum emission at -5cr, 5cr, 15cr, 25cr, and 60cr rms levels with the cr = 0.04 Jy beam“^. The yellow circle denotes the primary 
beam of SMA at 1.3 mm. Red circles denote the peaks of hot core (HC), mm2, mm3a, mm3b, the southern region (SR), north-east 
clump (NE), the NW and SE parts of the high-velocity outflow (OEIN, OEIS). Yellow crosses denote the BN object and the compact 
ridge (CR). The beam in the bottom left comer is from the SMA-only data. Colourmap shows the brightness (in counts per 400 s) of 
shocked H 2 emission ( Nissen et al.|2007| ). Corresponding to the yellow square, the insert panel shows the HC continuum from the 
SMA at 1.3 mm in detail as colourmap, with black contours of the 865 jum continuum emission at a spatial resolution of ~ 300 AU 
from Beuther et al. ( 2004[ ) overlaid. 


CDMS/JPL intensity is > 10“® nm^MHzl^ 


^ This is a coarse cutoff for identifying the strong emission lines, 
corresponding to spontaneous emission rate A^i ~ 10“^-^ Taking 
the frequency resolution of our observations into consideration (~ 0.8 
MHz), we found that 50% number possibility at certain rest frequency 
of a particular line has CDMS/JPL intensity > 10“^. 
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Table 1: Properties of substructures, which are denoted on the continuum image obtained from SMA observations at 1.3 mm 
(Figure]^. 


Source 

R.A.(a) Dec.(<5) Iv“ T„t N^,_i 

[J2000] [J2000] (Jy/beam) (K) (lO^^cm-^) (lO^^cm'^) 

5h35m~ -5‘>22'~ 3.49" x 2.57" 5.55" x 4.35" 

Alternative Designations 

HC 

155 ±16*1 23.7 ±4.6 

121 ±16*2 31.0±6.7 

14.6s 31.0" 2.75±0.3 126±13*3 29.6±5.8 11.9±0.7 

225 + 19*“ 16.3 + 2.9 

199±45*5 18.4±5.1 

CaA SMM3', CntB*^ 

hotcore^ 
SMAl^ 
Source F 

Source 

14.58s 30.97" 

14.46s 32.36" 

14.51s 30.53" 

14.35s 33.49" 

mm2 

108 ±4*1 11.0±1.7 

14.1s 27.1" 0.87±0.10 112±7*3 10.5±1.8 9.3±0.1 

174 ±14*“ 6.7 ±1.2 

173 + 40*5 6.7 ±1.9 

NWC*-', CcK SMM2', 

WC", CntD*^ 

cr' 

14.0s 36.9" 

SWC, SMMl' 

mm3 a 

103 + 3*1 12.8 ±1.6 

14.3s 37.9" 0.96 + 0.09 101 + 5*3 13.1 ±1.8 6.2±0.07 

159 ±12*“ 8.1 ±1.3 

101 ±26*5 13.1 ±3.7 

CbV, CntC^ 

mm3b 

88 ±6*1 15.9 ±2.7 

14.1s 36.4" 1.01 + 0.11 102 + 5*3 13.5 ±2.1 6.7±0.05 

160 ±17*“ 8.5 ±1.6 

89 ±24*5 15.6 ±4.8 

Cb2-'', CntC^ 

NE 

ITSs \9A" 0.38 ±0.09 43±3*i 20.1 ±4.2 3.3 ±0.06 

“CntI? 

SR 

14.0s 45.4" 0.41+0.06 43 + 3^ 14.2 + 2.9 1.0 + 0.02 

Cd^’, MM23"' 

OFIN 

13.4s 20.5" o- = 0.04 43 + 3^ <4.2 ~ 1.1 

PKr 

OFIS 

15.3s 42.5" o- = 0.03 43 + 3^ <3.1 -0.6 

PK2^ 

BN^ 

14.11s 22.7" 



Notes. Positions with upper case are from the references listed in the footnotes, without are from 1.3 mm SMA continuum observations. 
a. Peak specific intensity ly is obtained from each substructure of SMA-only continuum. 
bi. Excitation temperature is from CH 3 CN rotation map fitting in Figured 
b 2 . Excitation temperature is from CH^^CN rotation map fitting in Figur^ 

^ 3 . Excitation temperature is from HCOOCH 3 rotation map fitting in Figurej^ 
b 4 . Excitation temperature is from CH 3 OH rotation map fitting in Figure ^ 
bs. Excitation temperature is from ^"^S 02 rotation map fitting in Figure j9 

c. Excitation temperature is assumed the same as that in NE from CH 3 CN. 

d. H 2 column density is calculated from the SMA-only continuum at substructure peak using Eq.^^ performed in Section ^ 


4.2.2 


e. H 2 column density is calculated fr om the conversion of combined SMA-30m at temperatures derived from HCOOCH 3 , performed in Section 

R eferences. /. 


Menten & 


Blake et al. 


Be uther e t al. 
1 1995} ; 


2004 


4.2.2 


19961; 


Tang et al.^ 2010|: 


Favre et al.H 2011a|; 

pp2 


Murata et al. 


Zapata et al. 


m. Eisner et al. 


Wright et al, 


20TIT 


2008 


1992 


Beckwith et al 


1978 
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Fig. 3: Spectra taken towards selected positions (denoted by red circles in Figure]^ after imaging the whole data cube. The spectra^ 
resolution is smoothed to 1.2 kms“^ after the combination of SMA and 30 m data, and strong emission lines are labelled. A blow-up 
figure of the spectrum in the HC is Figure Al. 
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Using the ab ove method, the detection of strong lines are 
shown in Figure |A1| In addition, we can also detect multiple 
transitions of low-abundance COMs in Orion-KL, owing to the 
high sensitivity and recovering the full intensity for spatially 
extended molecular emission from the combined dataset. We 
use a synthetic spe ctrum fitting program ( |Sanchez-Monge|201 1[ 
Palau et ar]|2011[ ) in conjunction with molecular data from the 
JPL/CDMS catalogues, to simultaneously fit multiple COM 
lines under the assumption that all transitions are optically 
thin, are in local thermodynamic equilibrium (LTE), and have 
the same line widtlf^ Given these assumptions, the fitted 
flux densities are not always a perfect match to the observed 
spectra. However, in many cases we are a ble to robustly identify 
multiple transitions of a COM (Figure |A2| ). We discuss the 
properties of the identified COM species further in Section [L2.2[ 


With these two approaches, we identified more than 160 lines 
that we assigned to 20 molecules (including 25 isotopologues) 
over the entire band; we labelled the unblended strong emissions 
in Figureand list all the lines in Table |Al[ This table includes 
the rest frequency, quantum number, symmetry label, and upper 
state energy for each transition. Lines marked with “t” are from 
molecules containing less than six atoms as shown in Figure 
and those marked with are from molecules we designated 
as COMs and carbon chains, which are shown in Figure 
Our observations have a frequency resolution of 0.8125 MHz 
(~ 1.2 kms“^ at 1.3 mm). Thus, some observed lines with broad 
linewidths can be attributed to multiple transitions, when com¬ 
pared to the laboratory-measured rest frequencies recorded in the 
“Splatalogue” database. If potentially blended transitions for one 
line are from different isotopologues but have similarly strong 
CDMS/JPL intensities, we compared the spatial distribution of 
this line to those of the known isotopologues. Then, we either 


confirm a match or indicate all possibilities with in Table Al 


If there are multiple possible transitions from the same isotopo- 
logue at this rest frequency, they are also marked with 


3.2.2. Line imaging 


In general, all the identified species show different spatial 
distributions, velocity widths (FWHM AV), and velocity at 
line centre (Vpeak) in each substructure (see Figure I A3 1 on 


-1 


average, Vpeak ~ 6-8 kms ^ in the HC; Vpeak ~ 8-10 kms 
in mm2, mm3a, and mm3b; Vpeak ~ 10-12 kms“^ in NE). 
However, transitions from a single isotopologue have the same 
spatial distribution (except for CH 3 OH). In this paper, we have 
chosen one of the strongest unblended transitions from each 
isotopologue and produced their intensity distribution maps 
from the SMA-30m combined data. In avoiding contami¬ 
nation from other strong lines, we integrate the intens ity of 
each line through 40%-60% of their profile range in Figure 


The synthetic fitting program we use here is for line detection, but 
not for finding the best fit of column density or temperature to the ob¬ 
servations. Since the synthetic fitting is performed with optically thin 
assumption, minor lines of a particular species which cannot fit as well 
as the majority may be optically thick. 

The chosen velocity ranges are slightly larger than the FWHM of 
the lines, when accounting for the small variation in the peak velocity 
between different substructures. 
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Fig. 4: Intensity distribution maps of transitions from molecules containing <6 atoms, which are detected in the combined SMA- 
30 m data at 1.3 mm (corresponding quantum numbers are listed in Table |Al[ with “t” marks). The intensities have been derived 
by integrating the line emission over the velocity range shown in Figure A3 The filled beam in the bottom left comer is from 
SMA-30m data. White contours show the continuum from the SMA-only data (at -5cr, 5cr, 15cr, 25cr, and 60cr levels). Purple 
contours show ±lcr levels of the molecular emission (solid contours indicate positive fiux, and dashed contours indicate negative 
fiux). Green lines in the SO and CO maps sketch the outfiow directions. The white crosses denote the BN object and the CR, and the 
black crosses denote the positions of the hotcore, Source I, Source N, and SMAl, as in Figure All images have different colour 
scales (in Jybeam“^ kms“^), increasing from black to white, which are optimized to emphasize the features in the distribution of 
each molecules. 
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* Channel maps of more CH 3 OH lines are shown in Figure [A 6 | 


Fig. 5: Intensity distribution maps of transitions fro m C QMs and carbon chains detected in the combined SMA-30 m data at 1.3 mm 
(corresponding quantum numbers are list in Table |Al[ with marks). The intensities have been derived by integrating the line 

The filled beam in the bottom left corner is from SMA-30 m data. White 


A3 


emission over the velocity range shown in Figure 
contours show the continuum from the SMA-only data (at -5cr, 5cr, I5cr, 25cr, and 60cr levels). Purple contours show ±lcr levels of 
the molecular emission (solid contours indicate positive fiux, and dashed contours indicate negative fiux). The white crosses denote 
the BN object and the CR, and the black crosses denote the positions of the hotcore, Source I, Source N, and SMAl, as in Figure ^ 
All images have different colour scales (in Jybeam“^ kms“^), increasing from black to white, which are optimized to emphasize the 
features in the distribution of each molecules. 
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The combined data at 220/230 GHz covers a field of view of 
52", which is greater than the extent of the source (not including 
the outflows), so our maps exhibit negligible artefacts caused 
by residual missing flux. The high quality of our data allows 
us to map molecules, such as ^"^802, O^^CS, HNCO, H^^CO, 
^^CO, and CH 2 CO for the first time with both high resolution 
and sensitivity to spatial information on all scales. After sorting 
by different groups, we present the intensity integration map for 
each identified molecule in Figures]^ andCompared to previ¬ 
ous sub-mm and mm (especially interferometric) observations in 
Orion-KL (listed in the following paragraphs), the main features 
of our molecular distribution maps are: 

1. All the identified sulphur- (S-) bearing species in our 
dataset are diatomic/triatomic containing no hydrogen. 
Most of them have extended distributions (except for 
SO 2 and ^"^802), exhibiting a V-shaped morphology that 
extends from the emission peak at HC, through mm3a, 
mm3b to mm2. Their velocities at line centre also vary 
among substructures from ~ 6 kms“^ (HC) through 8-10 
kms“^ (mm2, mm3a, mm3b, and 8 R) to 11-12 kms“^ (NE). 



OC 8 and O^^C 8 (carbonyl sulphide) peak at 8MA1 
and have emission that extend to the CR. In addition, 
OC 8 is also detected in NE and 8 R. Although only one 
line of O^^C 8 is detected in our data (19 ^ 18, Eig,[^, 
other lines of this isotopologue have been previously 
confirmed and imaged by single-dish observations (e.g.. 


8 utton et n^|1985t Ziurys & McGonagle||1993[ [Tercero| 


et al.||2Q10 ). In addition to the detections by single-dish 

surveys (e.g., | 8 utton et al.||1985[ |Tumer |1991[ 8 chilke| 
et al.|1997 Tercero et al.|2010| )7previous interferometric 
surveys of OC 8 also reproduce the same V-shaped 
widespread intensity d istribution of this molecule, 
peaking at th e HC (e.g., [Wright et al.||1996( [Eriedel &] 
8nyder|2Q08| ). 


- 8 O 2 and ^"^802 (sulphur dioxide) peak at the hotcore 
and 8MA1 and exhibit extended emission in the NE 
direction, with Vpeak ~ 9-11 kms“^. Unlike the other 
8 -bearing species, no emission has been detected from 
these species in 8 R, and they do not display a clear V- 
shaped feature. However, there is a clear second velocity 


0-1 


component in the line profile of ^^802 (at ~ 22.7 kms 
Eigurej^ and 8 O 2 (at 22.1 kms“^ Eigurej^. In addi¬ 
tion, red- and blue-shifted ^802 trace outflows in the 
NW- 8 E direction (Eigure|A4|). 8 ingle-dish studies (e.g., 
80?; p^ke et al. | |1987i j Persson et al " 


et al.T2013 


■3T 


8 O 2 ; Esplugues et al.[2013|) conclude 


that this species is a good tracer not only for the warm 
dense gas, but also for regions affected by sho cks. While 


based on th e velocity around 22 kms ^ [Esplugues 


et al. 


( 2013F suggests that this velocity component 


may be associated with shocks caused by the explosive 
(NW- 8 E) outflow and/or the BN object, our data clearly 
show that this velocity component of 8 O 2 is spatially 
offset from BN and seems to be associa ted m uch more 
with the main NE- 8 W outflow (Eigure [A5[ ). On small 
scales, interferometric observations of Orion-KL also 
detect a compact 8 O 2 emission peak ( [Wright et al.|1996[ 
Zapata et al. 2011[ ) and strong extended emission of 


12 


V 


peak ■ 


20.5 kms ^ in Esplugues et al. {20131 


Fig. 6: Lines that only have one transition with low specific in¬ 
tensity (< 10 Jy beam“^) in our data, fitted with synthetic spec¬ 
trum from their emission peaks. 


34 


8 O 2 to mm 2 ( Beuther et al. 2005| . 


- 80 (sulphur monoxide) has multiple strong emission 
peaks towards HC, NE, mm3 a, and mm 3b an d exhibits 
strong self-absorption in mm2 (Eigure \A3\ . There is 
also significant emission extending along both the 
plateau of high- (Vpeak ~ 8-10 kms“^) and low-velocity 
(Vpeak ~ 7-8 kms“^) outflows. 8 imilar to ^"^ 802 , its line 
wings are broad (aV ~ 23 kms“^), and a second velocity 
component is c lear in the line profile at ~ 21 kms“^ 
(see Eigure I^A3I mm3b), which may be associated with 
shock event These morphological features (elongated 
emission peak(s) and broad line wi ngs of 80) have also 
been detected by both single-dish ( [Persson et al.[[2007[ 
Esplugues et al.poT^ and interferometric ( [Wright et al. 
1996[ ) line studies. 


- ^^C 8 (carbon monosulphide) has only one transition in 
our data set, but it is very strong, peaking at 8MA1, with 
a similar distribution to O^^C 8 -extended emission from 
NE through HC to the CR. Owing to its high dipole 
moment, th is molecule is consider ed to be a dense gas 
tracer (e.g., [Bronfman et al.[[1996|), and its detection in 
bot h CR and NE is con firmed by [Persson et ak] ( [20Q7[ ) 
and ffecero et al.] ( [ 2010 [ ). 


2. The distribution of simple CxHyOz (x,z=l,2...; y=0,l,2...; 
x-f y -f z < 5) molecules are co-spatially extended. 

- CO, C^^O, and ^^CO (carbon monoxide) have a common 
morphological peak at the HC (though the integrated 
intensity of CO peaks north of the HC d ue to its strong 
self-absorption in this source, see Eigure [A3[ ). Velocity 


The line wing of SO at ~ 7 km s~^ is broad, and its second velocity 
component is blended, so we c anno t tell the origin of the shock precisely 
from our observations, Figure 


A5 
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at line centre Vpeak ranges from 6 kms ^ (HC) through 
8 kms“^ (mm2, mm3 a, mm3b, and SR) to 10 kms“^ 
(NE, OFIN, and OFIS), with broad red- and blue-shifted 
line wings trace both high-velocity (CO , ^^CO) and low- 
velocity (C^^O) outflows (Figure A4). As the second 
most abundant species in the ISM, this species has been 
detected in numerous other observations of Orion-KF 
(C^^O;|Beuther & Nissen|2Q08l |Plume et al.B2012| CO; 


Zuckerman et al.|197^ Kwan & Scoville|1976|[Erickson 


et al. 1982[ Blake et al 


1987| (single-dish), Masson et al 


|Chernin & Wright|1996t |Zapata et al. 12009 [ Peng 

2012b (interferometry), ^^CO; Wilson et al. 2011 


1987 


et al. 


(single-dish)), showing the same morphology as we see: 
extended structure, with bipolar line wings covering 
from the NE through HC to the CR and tracing outflows. 

CH 2 CO (ketene) has a single transition ( 111,11 ^ lOijo) 
detected (Figure [^, peaking at mm3a and mm3b, with 
strong emission extending to the OFIS and the NE 


filament. With V, 


peak 


8 kms ^ and aV ~ 2.5 kms 


m mm2, mm3a, and mm3b, this line has also been 
confirmed by Sutton et al. ( 1985| ). Some ot her lines of 
this species are also detected in Orion-KF by [Lee & Cho 


( [2OO2I ). 

H^CO (formaldehyde) exhibits an elongated peak 
stretching from SMAl towards mm3 a. Similar to the 
simple organic CH 2 CO, it also shows extended structure 
along the SE lobe of the high-velocity outflow and NE. 
Single-dish detections (e.g., Persson et al.|2007| ) confirm 
that this molecule is present in both the CR and the 
outflows. 


HNCO (isocyanic acid) is the simplest organic molecule 
containing C, H, O, and N elements. Peaking at SMAl, 
it shows extended emission towards all the substructures, 
with Vpeak ranging from 6 kms“^ (HC) through 7-8 kms“^ 
(mm3a, mm3b, SR, and OFIS) to 9-10 kms“^ (mm2, NE, 
OFIN), and a second velocity component at ~ 19.8 kms“^ 
(Figure [7]). A lthough the formation pathway of HNCO is 
still not clear jHasegawa & Herbst|199'3|[Turner et al.|1999[ 


[Garrod et al.||2008t [Tideswell et al.||2010| ), its distribution 
is observ ed to be correlated with shock tracing mol ecules 
SiO (e.g., Zinchenko et al.||2000| ) and CH 3 OH (e.g.. 


& Turner|2005| ), suggesting that it is also produced in shocks. 


Meier 


In contrast, highly saturated COM^ are the typical COMs 
detected in HMSFRs, and they exhibit a variety of morpho¬ 
logical structures in Orion-KF. In general, N-bearing COMs 
have a similar peak to HNCO, at or near HC; while most of 
the 0-bearing COMs are co-spatial with CH 2 CO-peaking at 
mm3a and mm3b, and extending even to southern CR (e.g., 
HCOOCH 3 , CH 3 OCH 3 , and CH 3 CHO). This spatial segre¬ 
gation of the N- and 0-bearin g COM peaks has also been 
reported by other single-dish ( Blake et al.| 1987[ Persson 
|et al. 12007] ) and interferometric studies with high spatial res- 
olution (e.g. , |Wright et al.|19^ Beuther et aL|2005[ |Favre 


|et al.||2011at |Friedel & Widicus Weaver||2012| ). In addition, 

there are several other 0-bearing COMs identified in this 
study, which have unique morphologies. To understand the 


A COM that has a chain of carbon atoms linked together by single 
bonds and has rich hydrogen atoms filling the other bonding orbitals of 
the carbon atoms. 




-20 0 20 40 

Velocity (km/s) 




Fig. 7: Gaussian fitting to the line profile of ^^CO, C^^O, HNCO, 
SO 2 , and ^"^802 extracted from NE. Outflow signatures are seen 
in the broad line wings and the second velocity component in the 
HNCO, SO 2 , and ^^802 line profile. 


spatial differentiation of COMs on small scales (Figure [^, 
we compared the distributions of these molecules between 
our observations to other recent interferometric studies 
(Table[^. 

- CH 3 CN and CH^^CN (methyl cyanide) are two abun¬ 
dant N-bearing COMs. In our dataset, we detect the 
J = 12 ^ 11 ladder of both isotopologues around 220 
GHz. Most of the CH^^CN lines are blended, but 
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lines with K = 2,4,5,6 are clearly detected in the HC 
(Figure [^. In particular, K = 2 line of CH^^CN and 
CH 3 CN are clearly co-spatial (Figure |A3| ), peaking at 
SMAl (Vpeak ~ 5-7 kms“^) and extending to mm2, 
mm3 a, mm3b (Vpeak ~ 8-10 kms“^), with > 3 cr detec¬ 
tion in NE (Vpeak ~ 10-13 kms“^), exhibiting a V-shaped 
morphology. CH 3 CN is detected in SR and OFIN (S), 
but not for CH^^CN. This V-shaped morphology of 


Beuther et al. 2005 [Eriedel & Snyder 

20081 Guelin 

et al.|2008| Eriedel & Widicus Weaver|20^ 

L2t|Peng et al. 


2012a| ). Metha nol has previously been detected in the 
Orion outflow ( Wang et al.|2011| ), and furthermore, line 


8 - 1,8 ^ 7 o ,7 is known to be a potential Class I maser that 
traces shocks ( |Fish et aL]|2011| ). As a result, the SE tail 
can be reasonably explained as a feature of the outflow. 
In addition, it is also suggested that the SMAl-Source 


CH 3 Cr 

(2005); 

s has also been detected m e.g., [Beuther et al. 
i Zapata et al. ([2011 ); [Widicus Weaver & Eriedel 

(e.g., Menten & Reid| 1995 [ [Beuther & Nissen[[2008 ), 

( 2012 ): 

i Bell et al.|(2014[), as well as the NE extension 

which may explain CH 3 OH and H^^CO peaking in this 
area. 


warm and dense g as tracer (e.g., |Kalenskii et al.|2000[ 


|Araya et al.||2005| ), with gas phase and grain surface 
chemical models predicting that it is only detected 
in an environmen t with elevated temperatures (e.g., 
[Purcell et al.||2006| ). We use the CH 3 CN ladder and the 
CH^^CN (122 ^ 112 ) line to estimate the temperature 
of the s ubstru ctures where these lines are detected (see 
Section|4XT]). 

CH 2 CHCN (vinyl cyanide) exhibits strong compact 
emission north of the hotcore, with Vpeak ~ 5 kms“^ and 
a slight extension towards mm2. The northern peak is 
also coincident with the one in higher resolution maps of 
this species from SMA at 865yum (jBeuther et al.|2005| ), 
CAR MA at 1 mm (jEriedel & Snyder|2QQ8| ), and PdBI at 
3 mm ( jGuelin et al.|2008j )! 

CH 3 CH 2 CN (ethyl cyanide) emission also makes a 
V-shaped structure. Unlike the single-peaked structure 
of CH 2 CHCN, CH 3 CH 2 CN exhibits a dual-peaked 
structure, with a first emission peak at the hotcore 
(~ 6 kms“^), and a second peak at mm2 (~ 9 kms“^). 
Although there are two velocity components in both 
mm3a and mm3b (6 kms“\ 10 kms“^), there is no 
detection in southern CR, which is consistent with |Blake] 
et al. ( 1987j ). This observed morphology is consistent 


with other interferometric observations at high spatial 


resolution (Wright et a 

. 19961 [Beuther et al. [20051 

Eriedel & Snyder [2008 

Guelin et al. 

20081 [Widicusj 

Weaver & Friedel|20121| 

^eng et al.|2013 



CH 3 OH (methanol) has multiple transitions in our 
data (Eigure |A2| ). We obtained the channel maps of 
all the unblended strong CH 3 OH transitions in our 
dataset and found two types of spatial distributions 
(Eigure |A 6 | ). One type of line clearly shows extended 
emission from mm2 and SR (7.2-8.4 kms“^, e.g., lines 
of 80,8 ^ 7i,6E, 8 - 1,8 ^ ^ojE, and 3 - 2,2 ^ 4 -i, 4 E) 
towards the high-velocity outflow OEIN(S), exhibiting 
a similar SE tail to H^^CO at 7-10 kms“^ which has 
not been detected before; whereas the other type of 
line (e.g., lines of 154,ii ^ 163 ,i 3 E, 224,18 ^ 2 l 5 ,i 7 E, 

and 195,14 ^ 204 , 17 A-) does not show such a SE 

tail. However, both of them show multiple velocity- 
dependent emission peaks: towards SMAl at 5 kms“^ 
mm3a at 8 kms“^, and mm3b at 10 kms“^. The 
multiple-peak structure is different from the tradi¬ 
tional segregated dichotomous distribution of COMs 
in Orion-KL (i.e., the main emission peak is either at 
HC or mm3a(b)), though it has also been observed 
in other interferometric surveys ( j Wright et al. 1996[ 


HCOOCH 3 (methyl formate) has the largest number of 
detected lines in our data (Eigure |A2| ). There are two 
types of HCOOCH 3 in the whole observational band, 
torsionally excited lines (y = 1 ) and lines in ground 
state (y = 0). Both of them have identical distributions 
(Eigure that peak at both mm3 a and mm3b, with a 
clear detection in the southern CR, which is similar to 
CH 2 CQ. This morphology has also been reported by 


[Wright et al. 

( 1996[); [Beuther et al.j ([2005 ); Guelin et al. 

(|2008[); [Eried 

el & Snyderj (| 

2008[); [Eavre et al. (|2011a[); 

[Widicus Weaver & Eriedel 

[( 2012 |), indicating that this 


species arises m a cooler and less dense environment 
than the HC. 

CH 3 OCH 3 (dimethyl ether) peaks towards mm3 a and 
mm3b and extends to the southern CR, similar to 
HCOOCH 3 . Unlike the other COMs, its emission in 
mm2 is considerably stronger than towards the HC. 
This morphology has also been imaged with other 
interferometric studies (|Liu et al.||2002[ jBeuther et al. 


2005 


et al. 


Guelin et al.|2008y Eriedel & Snyder|2008t [Eavre 


2011bt|Brouillet et al.|2013|) 


CH 3 CH 2 OH (ethanol), the functional isomeif^ of 
CH 3 OCH 3 , exhibits extended emission from the vicinity 
of HC to mm2, mm3a, and mm3b, indicating the 
same distinctive V-shaped molecular emission as the 
N-bearing COMs. However, unlike CH 3 OCH 3 and the 
other COMs, it peaks in between the hotcore and mm3 a. 
The difference in the location of their peaks is exactly 
the same as reported by [Guelin et al.j (|200^; [Brouilletj 
et al.[ ( [2013 ) from PdBI and [Eriedel & Snyderj ( |2008[ ) 
with the EVLA and CARMA at high spatial resolution. 
In the grain mantle chemistry models, both isomers are 
produced from t wo radicals involved in the synthesis 
of CH 3 OH ice ( [Garrod et al.|[200^ . Therefore, their 
spatial differentiation may result from the chemical 
differences between the CH 2 OH and CH 3 O radicals 
across Orion-KL, e.g., variations in production rates or 
chemical evolution ( [Brouillet et al.|2013] ). 

CH 3 COCH 3 (acetone) is an 0-bearing COM, though, it 
has a similar distribution to CH 3 CH 2 CN, which peaks 
at SMAl, extends to mm2, and has no detection in the 
southern CR. Although its nature as a hot molecular 
core tracer has been confirmed by detections in both 
Orion-KL ([Eriedel et al.[[2()05 Eriedel & Snyder[[2008[ 


Widicus Weaver & Eriedel[[2Q12t [Peng et al.[[2Q13[ ) and 


Two compounds with the same molecular formula but different 
functional groups (radicals) are called functional isomers. 
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Sagittarius B2 (N-LHM) ( [Snyder et al.||2002| , whether 
this species mainly forms on dust grains or in the gas 
phase is still not clear. 


- CH 3 CHO (acetaldehyde) is tentatively detected, as all 
lines except one at 220.446 GHz (Table AT] ) are blended 
with HCOOCH 3 . Morphologically, this line has a similar 
distribution to CH 2 CO, peaking at mm3a and mm3b, 
with extended emission to the NE and S E tail. Although 
its origin is still unclear ( |Charnley|2004| ), it has b een de- 
tected in Orion-KL with single-dish telescopes ( jTurnerj 
|1989| ). Previous interferometric observations failed to 
detec t it because of filtering out the extended structure 
(e.g.,|Friedel & Snyder|20()8l jWidicus Weaver & Friedelj 

2012j ). However, its detection has been reported by 

recent CARMA and ALMA observations (Loomis, 
priv. comm), which shows the same distribution as our 
single-dish complemented interferometric observation 
result. 


- There are also several other organics, such as 
CH 2 OHCHO (glycolaldehyde), (CH 20 H )2 (ethy¬ 
lene glycol), and HCOOH (formic acid), which we do 
not detect in Orion-KL at 1.3 mm in the SMA-30m 
combined data. Of the first two, although there are some 
candidate detections from our synthetic fitting, they 
are blended with ot her nearby lines, and so cannot be 
confirmed. Although [Liu et al.j ( |20Q2| ) claimed to detect 
HCOOH in th e CR, in their data the line at 1 mm is only 
a 2 (T feature. [Widicus Weaver & Friedel| ( |2012| ) also 
did not detect any compact concentrations of HCOOH 
in our studied region either, and they conclude that the 
HCOOH emission must be extended on > 1700AU 
scales. 


- HC 3 N and HC 7 N are two cyanopolyynes (HC 2 n+iN, n is 
an integer). We confirmed the detection of HC 3 N(y 7 = 1) 
to the north of the hotcore and also tentatively detect 
HC 7 N. Although only two out of four of the HC 7 N 
lines predicted by the synthetic fitting programme 
are detected in our observation^^ these lines are 
unblended and peak at a similar location to HC 3 N. 
The excited-vibrational state {vq = \) of HC 3 N shows a 


similar confined structure north-east of the HC in Wright 
et al.j ( p^96| ); [Beuther et al.j ( |2005| ) and [Zapata et al. 


( [2QH[ ), where it is said to form a shell around the HC. 
HC 3 N is an especially useful dense gas tracer due to 
its relatively low optical depth tran sitions, large dipol e 
moment, and simple linear structure ( Bergin et al.|1996[ ). 
The HC 3 N(y 7 = 1) lines are believed to be primarily 
excited by a strong m id-infrared radiation field (e.g., 
[Goldsmith et al.[[1982 ). A gas-phase chemical model of 
a hot molecular core with high but transitory abundances 


of cya nopolyynes has been presented in [Chapman et al. 
( [2009[ ), where these species are formed from reactions 
between the precursor C 2 nH 2 (e.g., C 2 H 2 , acetylene) 
and N-bearing species via gas phase chemistry (e.g., 
C2nH2 + CN ^ HC2n+lN -f H). 


4. Rotational temperature and abundance estimates 

The above qualitative comparison of the morphology of differ¬ 
ent species show the chemical variations in the above-mentioned 
substructures. A more precise approach to studying these varia¬ 
tions is to measure the quantitative differentiations in chemical 
parameters towards individual substructures, i.e., the tempera¬ 
tures, molecular column densities, and abundances. 


In brief, from our detections, N- and 0-bearing COMs not 
only show two segregated peaks, but also exhibit more com¬ 
plicated spatial distributions. Most of the N-bearing species 
(except CH 2 CHCN) are co-spatial and peak at the hotcore. 
While some 0-bearing species (HCOOCH 3 , CH 3 OCH 3 , 
and CH 3 CHO) peak at mm3a and mm3b, and others present 
more unusual features (i.e., CH 3 OH has velocity dependent 
peaks at both hotcore and mm3a/mm3b; CH 3 CH 2 OH peaks 
in between mm3a and hotcore; and CH 3 COCH 3 peaks in a 
similar location to N-bearing COMs). This challenges the 
traditional hypothesis that all of these molecules are forming 
in the same chemical and physical conditions ([Garrod et al.[ 
[20081 [i:aaretal.| 2011 [ ). 


Unlike the saturated COMs, long carbon chains (organics 
with carbon-carbon multiple bon ds) are less like l y to be 
detected in star-forming regions ([Thorwirth[[2001[ [Rice &[ 
[Maier|2013[ ). 


- C 6 H is the parent radical of C 6 H (the first molecu¬ 
lar anion identified in the gas phase interstellar and 


circumstellar media; [McCarthy et al.[ 

2006), after its 

presence was predicted by[Herbst[([1981 

). It had not been 


detected in Orion-KL before. In our dataset, there is one 
unblended line at 230.109 GHz with J = 165/2 ^ 163/2, 
although most transitions identified by the synthetic 
fitting of C 6 H are blended. Emission from this tentative 
detected species appears to be confined to the HC. 


4.1. Temperature estimates 

Temperature is commonly used as a benchmark for comparing 
model results and observations. Since N-, 0-, and S-bearing 
species show different spatial distributions, we estimate their 
temperatures in each substructure separately. 


4.1.1. CH 3 CN 

A good tracer of temperature in high-density gas (n > 10^cm“^) 
is CH 3 CN, which is a symmetric-top molecule with no dipole 
moment perpendicular to the molecular axis, and its K ladders 
in rotational levels are excited solely by collisions. Provided that 
the collision rate is sufficient to thermalize the rotational levels 
within each K ladder, the rotational excitation temperature 
of CH 3 CN is equal to the gas kinetic temperature ( Boucher 
et al. 1980[ Wright et al.[[1995 ). Our observation band covers 
the J = 12 ^ 11 transitions of CH 3 CN around 220 GHz, with 
upper state energy levels from ~ 70 K to ~ 500 K, yielding a 
wide range for determining temperatures. If we assume that 
the emission from all of these transitions originate in the same 
parcels of gas along the line of sight in LTE, having the same 
size, velocity, and a single temperature, then rotation al tem- 
peratures can be derived from the rotational dia grams ([Hollis 


1982[[^ren & Mundy[[1984] [Qlmi et ar|[1993[ [Goldsmith & 
Langer|1999 ). Traditionally, in the optically thin limit, the level 


Lines with J = 196 ^ 195 and J = 205 ^ 204 are detected, but nei¬ 
ther J = 195 ^ 194 nor J = 204 ^ 203 are detected, so these potential 


detected lines are marked with in Table Al 
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Species 

Identification 

Peak 

Column density 
@Peak (cm“^) 

Other interferometric detections 

CH 3 CN, CH^^CN 

V 

SMAl 

1.4(17), 9.4(14) 

a, b, e, i 

CH 2 CHCN 

V 

North to hotcore 

>5.0(15) 

b, d, e 

CH 3 CH 2 CN 

V 

hotcore 

>1.1(16) 

a, b, d, e, i, k 

CH 3 OH 

V 

SMAl & mm3 a 

>2.5(17) 

a, b, d, e, f, i, m 

HCOOCH 3 

V 

mm3 a & mm3b 

>9.7(16) 

b, d, e, g, i, k 

CH 3 COCH 3 

V 

SMAl 

>1.5(16) 

c, i, k 

CH 3 OCH 3 

V 

mm3 a & mm3b 

>2.0(17) 

b, d, e, h, i, j 

CH 3 CH 2 OH 

V 

near SMAl & mm3 a 

> 1.5(16) 

d (x), e and j (only 1 line) 

(CH20H)2 

X 

- 

- 

- 

CH 20 HCHO 

X 

- 

- 

g. (some lines < 3cr) 

HCOOH 

X 



d (x), 1 (< 3cr), i (only from the ex¬ 
tended configuration) 

CH 3 CHO 

* 

mm3 a & mm3b 

>1.3(16) 

d (X), i (X), n (ALMA) 

CgH 

* 

SMAl 

>1.0(15) 

- 

HC 7 N 

* 

North to hotcore 

>1.0(16) 

- 


Wright et^ jl996}; b. B^her et al.||2QQ5' 


»nyder|j 2 UU^ e.|(juelin et al. 

" pOB ~ Liu et al. ' 2002} ; m. 


References, a. 


h. Favre et al 


201 lb ; i. Widicus Weaver & FriedelH2012 


c. Friedel et al. 12005); d. Friedel & Snyder 


j. Brouillet et al. 2013 ; k. Peng et al. 12013 


2008|; f. Z apata et al. 2011 ; g. Favre et al. 201 laj; 

Peng et al. 2012a| ; n. R. Loomis, private communication. 


Table 2: Identification of saturated COMs and complex carbon chains in Orion-KL. Confirmed COMs are marked with “V’’’ 
tentative detections with and nondetections with “x”. At each of their main emission peak (Figure [^, column density of 
CH 3 CN is obtained with optical depth correction, while column densities of the rest molecules are obtained at the temperature of an 
adjacent substructure from synthetic fittings, with an assumption that the species lines are optically thin (we denote them by 
x(y) = xxl0y. 


populations are then directly proportional to the line intensities 
of the K components. The parameters used to estimate the 
level populations of CH 3 CN transitions (i.e., the line strength 
Suij^^ (Debye^) and t he lo wer state energy El (cm“^)) in our 
data are listed in Table lA2l 

The K = 3, 6 lines are from ortho (o-) transitions, which 
have a statistical weight that is twice as strong as those of the 
other para (p-) lines. Radiative transitions between o-/p-types 
are usually prohibited, so we exclude these lines in the optically 
thin rotation diagram fitting (shown in black lines in Figure 
also see |Andersson et al.|1984| ). 

Panel 1 of Figure shows the synthesised fitting under the 
assumption of optically thin emission in LTE for the HC with 
the dashed lines. To reduce the number of free mrameters, 
we also assume that all the K = 0 - 8 component^ have the 
same line width. A satisfactory fit is obtained to all lines except 
K = 0, 1, 2. Since all transitions were observed simultaneously, 
calibration errors due to poor pointing or amplitude scaling 
are not responsible, and the remaining explanations for the 
deviations are as follows: (1). The K = 0, 1 lines are so close 
in frequency that they are blended in our fitting by the uniform 
line width we assumed. (2). It is also likely that lines at K = 2 
are more optically thick. 

Compared to CH 3 CN, lines of its rare isotopologue CH^^CN 
should be optically thin. However, CH^^CN (122 ^ II 2 ) is the 
only transition that is not blended with the broad line wings of 
CH 3 CN in substructures other than HC. Even in the HC, only 
the K = 2, 6 , and 7 lines are unblended (Panel 7 of Figure}^. 
Excluding the o-type lines (K=3, 6 ), the rotation temperatures 
in HC we derive from the main and rare isotopologues of 
CH 3 CN have large differences (black fittings in Panels 2 and 


The K = 9 line is blended by strong ^^CO (2^1) emission. 


8 of Figure [^. Optical depth seems to be the most likely 
reason for such inconsistency, so to investigate this hypothesis, 
we use the following iterative approach to correct the optical 
depth of CH 3 CN lines and derive the opacity-corrected rotation 
temperature Trot,p (p=HC, mm2, mm3a, mm3b, NE) of each 
substructure (red fittings in Panels 2-6 of Figure}^: 


(1) Obtain the rotation temperature Tl^ from weighted least 
squares fit to Boltzmann diagram, assuming that all the p-type 
CH 3 CN lines are optically thin (see fit to the black dots in 
Panels 2-6 of Figure}^; 

(2) Estimate the optical depth of CH 3 CN (122 ^ II 2 ) 
tk,i (K=2) by comparing its main beam temperature with that 
of CH^^CN (122 ^ II 2 ), and calculate the column density of 
line 122 ^ II 2 at Ti^p with the optical depth correction. Then, 
estimate the total column density of CH 3 CN as Nl^ assuming 
LTE (see Section [4^ for more details about the correction); 

(3) With Ti^p and as trial input parameters, estimate the 
optical depth for the remaining transitions in the CH^CN ladder 
^K,i (K=Q-8) with the RADEX code packag^^ (van der Tak 


et al. |2007| ), then correct their column densities line by line 


(dots and triangles in red in Figure}^ and obtain a new rotation 
temperature T 2 ,p; 

(4) Repeat step (2)-(3) until (\Tn,p - Tn-i,p\)/Tn-i,p < 10% 
(n=2,3,...; see fittings in red in Figure}^. 


RADEX is a one-dimensional non-LTE radiative transfer code, 
providing an alternative to the widely used rotation diagram method. 
Without requiring the observation of many optically thin emission lines, 
it can be used to roughly constrain the excitation temperature in addition 
to the column density. 
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Fig. 8 : Panel 1 ): CH3CN spectrum towards the HC in black, LTE model fit in dashed purple, and fitting with optical depth correction 
in filled red. Panels 2 - 6 ): rotation diagram of CH3CN derived for each substructure, with dots for P-CH3CN and triangles for o- 
CH3CN. Black fittings are only for the p-lines when assuming they are optically thin, and red fittings are for both o- and p- lines 
by optical depth correction. Panel 7 ): CH^^CN spectrum towards the HC in black and the LTE model fit in red. Panel 8:) rotation 
diagram of CH^^CN in HC with black dots (p-) and triangles (o-) from observation, and fitting is o-lines excluded by assuming lines 
are optically thin. Eor lines without detection, an upper limit equal to the 3cr rms is shown as arrow. The estimated temperature and 
total column density from each fitting is indicated in each panel. Errors are derived from scatter in the data points. 


The above approactj^ can be applied to sources HC, 
mm 2 mm 3 a, mm 3 b, and NE, where the CH3CN ladders and 
CH^^CN (122 ^ II2) can be clearly detectecj^ In SR and 
OEIN(S), the CH3CN lines are below the 3cr dkection thresh¬ 
old, so we assume the rotation temperature in these positions to 
be the same as that in NE. 


and mm 3 a(b), and tail towards the SE), which may be affected 
by their different temperatures. Since multiple transitions of 
CH3OH are available in our data, we can use them to derive a 
third temperature set. 


4.1.2. HCOOCH 3 , CH 3 OH, and 34 s 02 


HCOOCH3 has the largest number of lines detected in our 
dataset. Erom satisfactory synthetic fitting (espe cially for lines 
in the torsionally excited state v = 1 , Eigure A 2 ), we can 
reasonably assume that most of its lines are optically thin. Line 
emissions from this species are not contaminated by N-bearing 
species owing to their different spatial distributions, so we use 
them to derive a second set of temperatures. 


On the other hand, CH3OH and H^CO have an unusual 
spatial distribution (dual peaks /extended emissions at both HC 

In reality, molecules do not distribute homogeneously in the clouds, 
and they are far from being homogeneous in density and tempera¬ 
ture. The ultimate best-fit requires a sophisticated, fully 3Dnd physico- 
chemo-LRT modelling, which is beyond the scope of our paper. Our 
iterative approach is between realism and feasibility of the modelling, 
which fits both o-/p- CH 3 CN well. 

The detection of CHJ^CN (122 ^ II 2 ) in NE is 4o-. 


Eor S-bearing species, their extended emission indicates 
they may be sensitive to the ambient shock environment, so 
their temperature in each substructure may be different from 
the other species. However, in our data, some of the S-bearing 
species have only a single transition (O^^CS, ^^CS), and other 
lines may not be optically thin (SO2, SO, and OCS). Because 
^'^S02 is the only S-bearing rare isotopologue that has multiple 
transitions in our observations, we applied temperatures derived 
from transitions of ^'^S02 to the other S-bearing species. 


Using the lines list in Table A 3 we derived the rotation 
temperatures in each substructure from HCOOCH3, CH3OH, 
and ^'^S02 (Eigures J^-III) and list these temperature sets in 
bold face in Table ^ Temperatures derived from HCOOCH3 
lines in the ground and torsionally excited states separately have 
some differences, but they do not deviate significantly from 
temperatures derived from all of its lines. Therefore, we use the 
“mixed” temperatures in the end. 
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4.1.3. The Other organics 


In addition, since different organic molecules have specific 
spatial distributions, we also derive the rotation temperatures 
for the species that have multiple transitions with significantly 
different Eu/k e detected in our dataset (Figures |^V-VI, listed 
in Table |A3| ). In the condensations HC, mm2, mm3a, and 
mm3b, rotation temperatures of HNCO are similar to those 
derived from CH 3 CN within uncertainties; and rotation temper¬ 
atures derived from CH 3 OCH 3 and CH 3 CH 2 OH are close to 
HCOOCH 3 (y = 1) within uncertainties. Therefore, temperature 
sets from CH 3 CN and HCOOCH 3 are applicable to the rest N- 
and 0-bearing species (including simple molecules and COMs) 
in estimating their column densities. 


Sull^ 


3hc 


Eu 


^(Trot)Q(Trot)(e ^BTrot -e kBTrot) 


-1 




"kgSOOK _p kgSOOKl-l 


Debye^ 


(5) 

( 6 ) 


For the emission lines, we assume the Rayleigh-Jeans ap- 


proximation (^^^ 1) is valid, therefore Eqs.|3 

andp 

1 can be 

simplified as (Zeng et al.|2006 Eq. 4.26, also see 

Wilson et al. 


[ 200 ^ , 


Although temperatures derived from different species at the 
same substructure vary, HC is consistently the warmest substruc¬ 
ture, while mm2, mm3a, and mm3b are slightly cooler. In the 
outflow regions, there are no sufficiently strong lines for species 
other than CH 3 CN to derive the temperature, which may suggest 
these regions are cooler, exhibiting similar chemical properties. 

4.2. Column densities 

4.2.1. Molecular column densities 

For all the species we detected, we calculate their column 
densities (from both main and rare isotopologues) from the 
integrated spectral line intensities as below. 


Nu 3kB 


NT,a = 


ke (e^BTrot -1)^ 
he ^(Tj-ot) 


/ 

/ 


TydV 

(cm ^) 

(7) 

Tvdv 

(cm“^) 

( 8 ) 


Assuming that the observed emission in each substructure 
is homogeneous and fills the combined beam, the integration of 
measured main beam brightness temperature within the velocity 
range f TB(r')dL' can be substituted for the last term in the above 
equation: 

T,ot J T,dv=j^^ JTB(v)dv (Kcms-1) (9) 


The optical depth (in velocity) along the line of sight is 
( |Zeng et al.|2006| Eq. A.25): 

Tv = ^AuiNu(eTO - 1 ) 0 ( 1 ;). (1) 


where c is the speed of light, Fb the Boltzmann constant, h the 
Planck constant, v the line rest frequency, 0(i;) the distribution 
function of the line shape in terms of velocity v. Here, Agi is the 
average spontaneous emission rate from the upper state Eg into 
the lower state Ei, and it is calculated from line strength Sui and 
dipole moment ji as 


3hc3 gu 


( 2 ) 


After integrating the optical depth over the observed 
linewidth, f (X>(v)dv = 1 , and the column density Ny of the line 
at y Hz is 


N 


U — 


87ry^ 1 



(cm ^) 


(3) 


When assuming LTE, the total column density of all transi¬ 
tions for a molecule a is 


NT,a = —Q(Trot)e‘'BTro. (cm-2) (4) 

Su 


where Q(Ti-ot) is the partition function for the given (rotational) 
excitation temperature of each sources, which is interpolated 
from a table of partition functions for fixed temperatures, 
obtained from CDMS/JPL. 


where Ta,o is the optical depth at line centre, and f T^{v)&v is 
measured from Gaussian/hyperfine structure (HES) fitting by 
using Gildas software package. 


The column densities and related uncertainties of species 
listed in Table || (CH 3 CN (CH^^CN), HNCO, HCOOCH 3 , 
CH 3 OH, CH 3 OCH 3 , CH 3 CH 2 OH, and 3^802) can be estimated 
directly from their rotation diagrams and the scatterings of the 
data points (Figure |9]). For the species whose transitions are 
not sufficient to derive the rotation diagrams, we estimate their 
column densities by assuming that transitions of related species 
have the same Tj-ot^ 

• N-bearing species have similar temperatures to those derived 
from the opacity-corrected CH 3 CN, 

• S-bearing species have similar temperatures to ^^S 02 , 

• 0-bearing species have similar temperatures to HCOOCH 3 , 

• Given the similarly unusual spatial distribution between 
H^^co and CH 3 OH, we also estimate the column densities of 
H^CO at temperatures derived from CH 3 OH. 


In general, we assume that the observed transitions are opti¬ 
cally thin {jafi < 1, ~ 1). However, for several molecules 

with large abundances, we observe the corresponding transitions 
of their main and rare isotopologues, e.g., OCS (J = 19 ^ 18)- 
Qi^CS (19 ^ 18), (2 ^ l)-CO (2 ^ l)-C^^O (2 ^ 1) and 

CH^^CN (122 ^ Il 2 )-CH 3 CN (122 ^ II 2 ). By measuring the 
ratio between main beam brightness temperature of the main 
line Tb, af,o and its rare isotopologue Tb,^,o, we can estimate 


the optical depth at line centre of a given transition Tq, 0 (Myers 
let al.|1983| ): 


Her e, S uifJ^ can be calc ulated from CDMS/JPL line intensity 
^(Trot) ( |Pickett et al.|1998| ). 


l-exp(-TQ,,Q/^^) ^ Tb, ^,0 

l-exp(-T^,o) ^ TB,a ,0 


( 10 ) 
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Fig. 9: Rotation diagrams of the N-, 0-, and S-bearing species with multiple unblended transitions detection in HC, mm2, mm3a, 
and mm3b. Panels in Row I show two states of HCOOCH 3 lines, with the blue and black dots from observation of ground state 
(y = 0) and torsionally excited state (v = 1) lines, respectively. With the optically thin assumption, fittings for each state are plotted 
in blue or black, respectively, and fittings for all the transitions are plotted in red. Non-detection of HCOOCH 3 lines are not shown. 
Panels in Row II-VI show fittings with the optically thin assumption, with 3 (T upper limits from non-detections plotted as arrows. 
Blue dashed lines in Row III are fittings with upper limits included. The estimated temperature and total column density from each 
fitting is indicated in each panel. Errors are derived from scatters in the data points. 
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Species 

HC (K) 

mm2 (K) 

mm3 a (K) 

mm3b (K) 

NE (K) 

CH 3 CN 

155 + 16 

108 + 4 

103 + 3 

88 + 6 

43 + 3 

CRi^CN 

121 ± 16 

- 

- 

- 

- 

HNCO 

185 + 52 

121 + 18 

105 + 25 

97 + 19 

- 

HCOOCH 3 

126 + 13 

112 + 7 

101 + 5 

102 + 5 

- 

CH 3 OH 

225 + 19 

174 + 14 

159 + 12 

160 + 17 

- 

CH 3 OCH 3 

71+7 

87 + 2 

102 + 6 

92 + 8 

- 

CH 3 CH 2 OH 

137 + 33 

73 + 15 

70+17 

68 + 11 

- 

^^*802 

199 + 45 

173 + 40 

101 + 26 

89 + 24 

- 


Table 3: Rotation temperature sets for each substructure derived from species, each of which has > 4 transitions with significantly 
different Eu/ke detected in our data. The substructures where there are not enough sufficient lines to derive the rotation temperatures 
are denoted by Temperatures in bold face are used to estimate the column densities and abundances of species whose rotation 
diagrams are unable to be derived. 


where is the intrinsic abundance of the main isotope (e.g., 
compared to its rare isotope (e.g., in the ISM (Table 
1^ e.g., [Wilson & Rood|1994t|Chin et al.|1996| ). 


Using the above method, we calculate the column densities 
of 25 isotopologues (for both optically thin an d th i ck ca ses) 
in each substructure, and present them in Tables |A5^ , |A 6 ^ and 
|A7t . In Table [^ we list the transitions for which we corrected 
their optical depths. We found that CO (2 ^ 1), OCS (19 ^ 18), 
and CH 3 CN (122 ^ II 2 ) are significantly optically thick 
(r = 6 - 117) in all substructures where they are detected. The 
validity of the assumption that the re maini ng molecular lines 
are optically thin is discussed in Section [4^ 


In addition, the main emission peak(s) of some COMs 
(CH 3 COCH 3 , CH 3 CH 2 OH, CH 2 CHCN, C 6 H, and HC 7 N) do 
not strictly coincide with either of the continuum peaks. In de¬ 
termining their column densities at their main emission peaks, 
we use their strongest transitions (shown in Figure and the 
temperature from an adjacent substructure. Lines are assumed to 
be optically thin, and molecular column densities could be un¬ 
derestimated if this assumption is incorrect, so the lower limits 
are given in Table 


4.2.2. H 2 column densities 


1987 ), extrapolated from jOssenkopf & Henning' { 1994] ); yu is 
the mean molecular weight of the ISM, which is assumed to be 
2.33; and mn is the mass of an hydrogen atom (1.67 x 10“^"^ g). 
The results are reported in column Nh 2 ,i ofTable[^ 

From the SMA-only observations, the total continuum flux 
at 1.3 mm is 22.6 Jy. H owever, extrapolating the f lux Sy from 
existing SCUBA data ( jPi Francesco et ah |2008| ) at 450 jim 
(28 700 Jy) and 850 jim (2 090 Jy), the total flux should be closer 


to ~ 400 Jy at 1.3 mm (assuming Sy oc Johnstone & Bally 


|1999| ). The difference implies that > 90% of the flux has been 
filtered out by the SMA observations. In contrast, (2 ^ 1) 
from a combination of 30 m and SMA observations is detected 
towards all the substructures, and it should not suffer from the 
same spatial filtering as the SMA-only continuum data. As a re¬ 
sult, it is a more reliable probe of the total column density of gas 
in the cloud at all spatial scales. Moreover, as a rare isotopologue 
line, it can be reasonably assumed to be optically thin, so we use 
it to make a second independent estimation of the H 2 column 
density in all the substructures. By assuming a constant abun- 


dance in all substructures 


__ Nh 2 ~ Nco Nh2 

( [Wilson & Rood||1994| ), H 2 column densities can be calculated 
from the column densities of optically thin emission (re¬ 
ported in the column Nh 2,2 of Table [^. 


In Orion-KL, at 1.3 mm wavelength, the continuum emission is 
dominated by dust-neutral gas emission ( [Beuther et al.[[2004| ). 
Assuming the dust emission at mm wavelengths is optically thin 
and coupled with the g as, which has a sin gle excitation tempera¬ 
ture, Tex ~ Trot = Tdust ( Hildebrand|1983[ ), we use the SMA-only 
data to calculate t he H 2 column densitie s for different substruc¬ 
tures, according to Schuller et al. ( 2009[ ) : 


Nho 


IvR 


Fy(T(iust) Qa 


(cm ^) 


( 11 ) 


where ly is the peak specific intensity for each substructure on 
the continuum (in Jy beam“^ taken to be the 3cr upper limit in 
the outflow regions); R is the isothermal gas to dust mass ratio 
(taken to be 150, from [Drame|[2011[ ); By(Tdust) is the Planck 
function for a dust temperature T^ust; is the solid angle of 
the SMA beam (in rad^); /Cy = 1.0 cm^g“^ is the dust absorption 
coefficient at 1.3 mm (aspming a model of agglom erated grains 
with thin ice mantlej^ for densities 10^ cm“^ ( Blake et al. 


At around 100-150 K, water ice starts to be evaporated. However, 
in the line of sight, the envelope of the protostar may not have the same 


4.3. Molecular abundances 

Subsequently, to study the chemical properties on the small 
scale of Orion-KL, we derived the molecular abundances with 
respect to H 2 (c onver ted from C^^O) in each substructure 
(Tables [A5tl,[A6tl, and[A7tl). 

Figure presents the molecular column densitiy and the 
abundance of each species with respect to Nh 2 , 2 - Generally, 
all the molecules have higher column densities in the central 
condensations (HC, mm2, mm3a, and mm3b) than the extended 
regions in the outflows (NE, SR, OFIN, and OFIS). While the 
abundances of some species have relatively small variations 
among the substructures (e.g., SO, CO, C^^O, and ^^CO), 
some vary substantially. Specifically, N-bearing species are 

temperature at different radii, so ice may not be completely evaporated. 
Therefore, without knowing the real geometry of the source, we as¬ 
sume the dust absorption coefficient according to the “thin ice mantle” 
ca se. Moreover, this as sumption is consistent with the modelling result 
in [Gemer et al.[p014| ). In the extremely case, if we assume the ice is 
completely evaporated, the H 2 column density derived from continuum 
will be 2-4 times lower than our current estimation. 
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significantly underabundant in the SR and the outfiow regions; 
0 -bearing organic molecules containing < 6 atoms (e.g., 
H^^CO, CH 2 CO, CH 3 OH) have more extended emission across 
all the substructures, with nearly constant large abundances in 
all the substructures except for NE. The S-bearing molecules 
are an intermediate case: most of them are not detected in the 
high-velocity outfiow, but are clearly detected in NE (OCS, 
SO 2 , ^"^802, ^^CS) and (or) the SR (SO), with lower abundance. 
Such abundance differentiation among species is consistent with 
the spatial differentiation shown in Eigurej^ 


4.4. Error budget 

In this section, we consider the effect that the various assump¬ 
tions that we made in our quantitative analysis have on the re¬ 
sults presented in this paper. 

1. Trot (rotation temperature): The traditional rotation diagram 
method is subject to a number of uncertainties owing to 
various assumptions, e.g, line emission is assumed to be 
o ptically thin and LTE is assumed to apply (see di scussion 
in |Snyder et al.|2005}|Herbst & van Dishoeck|20Q9| ). If these 
assumptions are valid, the points on the population diagram 
will be well-fit by a straight line, with the slope proportional 
to the negative reciprocal of the temperature. However, 
this is not the case (see Eigure [^. According to Table [3 
the volumn density in each substructure is ~ 10 ^ cm^ 
(assuming a source size ~ 10 times of the beam size, see 
also |Blal^ et SI]|1987t [Plume et"ST]|2Q12| ), so LTE can be 
a reasonable assumption (except for such extreme cases as 
masing and sub-/super thermal excited conditions), and the 
main reason for the fittings with uncertainties is the optical 
depth of the used lines. 

(1) In Section Em we obtain the optical depth of 
CH 3 CN (122 ^ 112 ) by comparing its observed main beam 
brightness temperature to that of CH^^CN (122 ^ ll 2 )- 
Then, we use this optical depth as a standard in non-LTE 
radiative transfer fitting and derive a lower Trot in each 
substructure. In the HC, the rotational temperatures derived 
from CH 3 CN (155 + 16 K) and CHf CN (121 + 16 K) are 
slightly different, but they agree within the uncertainties, 
and column densities of N-bearing species are insensitive 
to temperature uncertainties of this magnitude (Table [A 6 ^ ). 
On the other hand, temperatures derived from HNCO agree 
with the above Trot within the uncertainties, so our approach 
should be applicable. 

(2) In Section Em T'rot for most O bearing species are 
assumed to be similar to those derived from HCOOCH 3 . 
In general, the temperatures derived from two different 
states of HCOOCH 3 lines are consistent, while the tor- 
sionally excited HCOOCH 3 (v = 1) lines appear to trace 
lower temperatures than lines in the ground state (y = 0 ), 
because LTE assumption may not be applied in sub-thermal 
conditions (see also jPemyk et al.||200^ . In addition, we 
derived a third temperature set from CH 3 OH, which is 
in general higher than those derived from HCOOCH 3 , 
CH 3 OC H 3 , and CH 3 CH 2 OH. Synthetic spectrum fitting 
(Eigure |A2| ) indicates that some lines of CH 3 OH may be 
optically thick, which may explain the higher measured Trot- 
On the other hand, CH 3 OH lines have broader line widths 
(~ 5 - 7 kms“^) than HCOOCH 3 lines (~ 3 - 5 kms“^) in 


the central condensations, implying that they trace gas with 
different temperatures. 

(3) Eor S-bearing species. Trot derived from ^"^802 varies 
from HC through mm2 to mm3 a and mm3b. However, 
detected lines in mm3 a and mm3b are not sufficient to 
derive a reliable temperature, because we lose the constraint 
of fittings from the non-detection line with high Eu/ke at 
229.987 GHz (Eigure |^II). Nevertheless, we estimate the 
temperature upper limits for both substructures by including 
the non-detected lines (dashed blue fittings), and find the 
differences for column densities derived at temperatures 
from 89-101 K (non-detection are excluded in fittin gs) t o 
176-193 K (upper limits, columns in italic font in Table |A7P 
are by a factor of < 3. The insensitivity of S-bearing column 
densities to this temperature regime may also explain the 
slightly higher T^t with large uncertainties in HC and mm2. 

2. Nt (molecular column density): In Section |4.2| we first 
assume that all the identified molecular lines are optically 
thin. Then, we use Eq. [T^ to obtain the optical depth of 
those main isotopologue lines whose rare isotopologue lines 
with the same transition are detected. Einally, using Eq. 
we correct the column density of the main isotopologue. 
In this approach, by assuming the rare isotopologue lines 
are optically thin, we list the optical depths and abundance 
ratios for (2 ^ 1), CO (2 ^ 1), OCS (19 ^ 18), and 
CH 3 CN (122 ^ 112 ) with respect to the corresponding rare 
isotopologue lines with and without optical depth correction 
in Table [U 

However, Eq.|^is based on the assumption that the main and 
rare isotopologue lines have the same excitation tempera¬ 
ture, which is often not the case. Therefore, a more accurate 
calculation requires the large velocity gradient (LVG) 
approach and proper chemical modelling. Nevertheless, 
our analytic calculations suggest that, the optical depth 
is the main situation driving the error in the estimated 
molecular column density. Eor abundant molecules such as 
CH 3 OH, HNCO, SO, and SO 2 , where we do not observe 
a corresponding transition of their rare isotopologues, our 
assumption that their lines are optically thin may result 
in underestimating their column densities by a factor of 
7-10 (coming from the average differences for OCS, CO, 
and CH 3 CN column density values before and after optical 
depth correction in Eigure [T^. 

Aside from the optical depth, the errors also come from 
our assumption of LTE. Assuming that the transitions of 
one species have ideal Boltzmann distributions at the same 
Trot, molecular total column density can be derived from 
one transition. However, as represented by points on their 
rotation diagrams, multiple transitions of certain species 
(e.g., 34 s 02 , HNCO, and CH 3 CH 2 OH, Figure||( show wide 
scatter, indicating that uncertainties in the estimation of 
Trot lead to the over-/under-estimation of the tot al column 
den sitie s by a factor of 0.5-2 (bold face in Tables [A5^ , j A 6 ^ , 
and |A7p . Eor the species whose total column densities are 
not able to be estimated from rotation diagram fittings owing 
to the lack of sufficient transitions, our LTE assumption 
(using Eq. may also result in the under/over-estimation of 
their total column densities of the same magnitude. 
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* Values of are derived from the temperatures the same as HCOOCH3; Values of are derived from the temperatures the same as CH3OH. 


* CeH is not given because of large uncertainties of its tentative detection. 




Fig. 10: Molecular column densities (panels in the left column) and abundances (panels in the right column) towards 8 positions in Orion-KL. The molecular 
abundances of each species with respect to H2 are converte d fro m C^^O, a ssuming a constant N(-i 8 o/Nh 2 ~ 1.79 x 10“^ ratio. Column densities and abundances 
with the uncertainties are measured as mentioned in Section [44] and Tables |A^A7| Coloured solid lines show the differentiation for each isotopologue calculated 
from optically thin assumption (arrows highlight the upper limits), while dashed lines show the differentiation calculated with optical depth correction from Eq.fToj 
For abundant molecules (e.g., HNCO, SO, SO2, CH3OH) which we have not measured their line optical depths, their densities could be underestimated by a factor 
of 7-10. 
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Species 



HC 

mm 2 

mm3 a 

mm3b 

NE 

S 

OFIN 

OFIS 

rco ( 2 ^ 1 ) 



54 

53 

117 

107 

74 

90 

36 

22 

NcOr/Nci8o 

t 

560 + 25 

1263.0±21.8 

1541.0+57.1 

2422.0+46.6 

2765.0+72.3 

3100.0+52.9 

2934.0+18.8 

2640.0+33.7 

1517.0+29.8 

Nco/Nc18o 

★ 


23.4+0.4 

29.1+1.1 

20.7+0.4 

25.8+0.7 

41.9+0.7 

32.6+0.2 

73.3+0.9 

68.9+1.4 

'^l^co ( 2 ^ 1 ) 



l.l 

0.25 

2 

1.5 

2 

1.5 

0.02 

0.08 


t 

7.2 ±0.9 

10 . 8 + 0.2 

7.3+0.2 

14.6+0.2 

13.1+0.2 

12 . 8 + 0.0 

IO. 6 + 0.0 

7.6+0.0 

9.7±o.3 

Ni3(2o/Nc18o 

★ 


6 .5+0.1 

6.5+0.2 

6 .3+0.1 

6 . 8 + 0.1 

5.5+0.0 

5.5+0.0 

7.5+0.0 

9.3+0.3 

rocs (19^18) 



5.5 

8.5 

6 

10 

<6 

<6 

<2 

— 

Nocs^/Nqiscs 

t 

77 + 7 

79.7+0.7 

81.1+2.0 

80.1+3.0 

78.2+0.0 

< 195.2 

< 121.3 

< 101.1 

— 

NoCs/NqIBcs 

★ 


14.4+0.1 

9.5+0.2 

13.3+0.5 

7. 8 + 0.0 

>32.4 

> 20.2 

>43.7 

— 

’"CHjCN ( 122 ^ 112 ) 



18 

11 

8.5 

9 

17 

<3.5 

<13 

<4 

NchjCNt/Nchiscn 

t 

77 + 7 

79 . 6 + 3.0 (76.6+2.4)* 

77.7+3.1 

78.1+1.9 

79.8+0.3 

82.3+6.7 

< 168.3 

< 176.3 

< 130.0 

NcHjCn/NchBcn 

★ 


4.4+0.2 (4. 3 + 0 . 1 )* 

7 .U 0.3 

9. 2 + 0.2 

8.9+0.0 

4. 8 + 0.4 

>46.6 

> 13.6 

>31.9 


* Ratio at temperature derived from CHj^CN. 


Table 4: Optical depth and abundance ratios for CO, ^^CO, OCS, and CH 3 CN with respect to their rare isotopologues in each 
substructure of Orion KL. The values are derived from the single strongest transition of each species using Eqs. and ^ "t or 
★ mark the estimations with or without optical depth correction. “> (<)” come from the 3cr limit of the non-detection of O^^CS 
(OCS) and CH^^CN . T he ratios between isotopes of ^^C/^^C and ^^O/^^O are provided from the local ISM ratio measured in 


Wilson & Rood (1994). Uncertainties on the measured values are typically < 10% (written as the subscript), as determined from 


Trot, partition function Q(Trot), and f TB(v)dv of the transition. 


In addition, f TB(v)dv measured from Gaussi an/H FS fittings 
brings on average < 10% uncertainties (Table [A4| ), which is 
negligible compared to the chemical differentiation amongst 
the subs truct ures (written as the subscript in Tables |A5^ , 
[A 6 t , and|A7t). 


3. Nh 2 and molecular abundance: Because no continuum 
information is available from the combined data, we used 
two methods to calculate the column density of H 2 (Table[^: 
from SMA-only continuum (Nh 2 ,i) and from SMA-30m 
C^^O (Nh 2 , 2 )- The molecular abundances were derived 
with respect to H 2 , 2 - Comparing the column density of H 2 ,i 
and H 2 , 2 , we found Nh 2 , 2 /Nh 2 ,i is ~ 0.9 in mm2, ~ 0.5 in 
mm3a and mm3b, ~ 0.4 in HC 0 ~ 0.2 in NE, but 0.1 in 
SR. Several reasons may lead to such difference: (1) the 
specific intensity (Jy beam“^) used to calculate the column 
density depends on the beam. The combined beam is larger 
than the SMA-only beam by a factor of ~ 2.5, so it is 
reasonable that for larger bea m size, we obtain low er H 2 
column density (similar case in |Beuther et al.|2007a ). (2) In 
producing the continuum, we select the “line-free” part, but 
there are still lines with low emission (Figure |A2| ), which 
may contaminate the continuum; (3) we simply assume 
that the the gas-to-dust ratio and the volume densities in 
all substructures are similar within an order of magnitude 
(~ 150, ~ 10^ cm“^), hence also the dust opacity. Since we 
do not know the size of each substructure, gas is possibly 
much denser in a warmer source (HC, mm2, mm3a, and 
mm3b) than in the other cooler substructures (NE, SR), so 
the homogeneous volume density assumption may not be 
accurate in the whole region. (4) Grain growth in the en¬ 
velopes can also be a reason for difference in dust e missivity 
from centre to extended substructures at (sub-)mm ( |Miotello| 


et af 


2014|). (5) We also note that column density 


in the central region is a factor of 30 higher than in the 
outflow. Judging that Nh 2,2 is derived with the assumption 
that N( 218 o/Nh 2 = const, in all the studied substructures, the 
variation of N(-i 8 o/Nh 2 ratio in different substructures can 


be another factor that leads to error in our estimation. 

4. In addition, another simplification in our calculation is 
an assumed constant beam-filing factor. Determination of 
molecular abundances in Orion-KL substructures relies 
upon an accurate source structure model for each molecule. 
When the source is smaller than the beam size, the main 
beam brightness temperature is diluted, leading to an 
underestimated brightness temperature for more compact 
structures. To correct this, an appropriate correction factor 
needs to be applied. Since we do not know the accurate sizes 
of the substructures, we simply assume that the filling factor 
is unity in all the calculations. 


5. Discussion 

These observations comprise a large set of molecular lines 
(over 160 lines from 20 species, including 25 isotopologues) 
from combined SMA and IRAM 30 m data towards Orion-KL. 
Resolving Orion-KL into eight substructures and comparing 
S-, N-, and 0-bearing molecules, we find complicated spatial 
morphologies for COMs and strong chemical differentiation 
over the observed substructures. 


5.1. Comparison with the other results 

Previous quantitative studies of Orion-KL mainly focused on the 
HC, CR (with mm3 a and mm3b treated as a single substructure), 
and plateau (OFIN, OFIS), therefore comparison of our esti¬ 
mated molecular abundances with previous results is limited to 
these substructures. 

1. Since different molecules have different chemistry, the 
rotation temperatures derived using different tracers in the 
same substructure vary. Comparing with the results of others 
studies, the variation in Tj-ot for different species can be 
explained. 


In the extreme case, the ratio is 0 . 8 - 1 , if assuming the ice mantle - Our opacity-corrected temperature in HC derived from 

is completely evaporated. CH 3 CN is 120-160 K, which is a factor of 1.5 lower 
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than reported by Goddi et al. ( |2011a| ) from NH 3 at a 
higher angular resolution (~ 0.8"). Given that HNCO in 
our dataset probes the similar temperatures, our slightly 
lower temperatures for N-bearing molecules may trace 
the envelope region of the HC. 

Our adopted temperatures fro m HCQQCH 3 agree with 
the temperatures estimated by [Favre et ^ ( |2011a| ) from 
observations of this species at similar spatial resolu- 
tio However, they are lower than those reported in 
Brouillet et al.| ( MT3] ) owing to our spectral resolution. 
Multiplets of CH 3 OCH 3 and CH 3 CH 2 OH are not 
resolved in this paper, so we treat each detected line of 
these species as a single component, and temperatures 
obtained from rotation diagrams are potentially underes- 
timatec 0 


Temperatures probed by CH 3 OH (> 200 K in HC, ~ 170 
in mm2, and ~ 160 K in CR) are higher than those 
probed by the other molecules, probably because of the 
lack of optical depth corrections, or because they trace 
gas with different temperatures (discussed in Section 
|4.4| part 1 (2)). Nevertheless, these temperatures agree 
with those derived from th e same species by |Sutton 


et al.|dl995| ) (^^CH 30 H) and |Friedel & Widicus Weaver 
( |2012| in HC and CR, and by Beuther et al. ( 2005| ) in 
HC and mm2. 


- We note that S-bearing species may be affected by 
shocks in the environment. Although they were not 
often derived previously, and the high Eu/ke transition 
leads to large uncertainties in the rotation temperatures 
derived from ^"^ 802 , our estimated tempe ratures are 
within the uncertainty of those derived by |Esplugues| 
leT^ ([ 20 T 3 I) from SO 2 in HC (190 + 60 K) and CR 
(80 + 30 K). 


2. Due to the compactness of the detected substructures, col¬ 
umn densities of CO, ^^CO, and C^^O are ten times higher 
than those from lower resolution obs ervations in single-dish 
observations of [Wlson et al.| ( |2011| ) (EWHP ~ 13") and in 


( 2012| ) (Herschel/HIEI beam ranges from 13"to 


Plume et al. 

40"). If taking the beam ratio into account, our results are 
consistent with previous measurement. 


3. In both HC and CR, the abundances of SO, SO 2 , HNCO, 
HC 3 N, CH 3 OH, CH 2 CH CN, CH 3 CH 2 CN, and CO are 
compa rable to the values in Blake et al. ( 1987| ); Wright et'af] 
( |1996| ), and |Esplugues et al.| ( |2013| )7within our factor of 


1.5-2 uncertainties. HJ,^C0 is lower than the value derived 


Blake et ak] ( |1987| ), assuming a standard 


from H 2 CO 

^^C/^^C ratio, but is comparable to the value in [Wright 


let al 


and C 


( 1996 ). After the optical depth corrections, our OCS 
i 3 CN abundances are ten times higher than those 


in |Blake et~aL] ( |1987| ), because the latter estimation was 
made assuming lines of these species to be optically thin. 


Favre et al. 


(2011 at , 


MFl-5 are the HCOOCH 3 emission peaks in| __ ___ 

which coincide with mm3b (MFl), HC (MF2), mm3 a (MF3), and mm2 
(MF 4-5) in our con t inuum . 
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Brouillet et al. 


(20131 resolves the CH 3 OCH 3 multiplets and uses 


them to estimate temperatures as well as column densities. However, if 
each multiplet is treated as a single component, it would derive temper¬ 
atures consistent with those we estimate. 


However, o ur abundances of PCS and ^^CS are comparable 
to those in |Tercero et al.| ( |2010| ), within our factor of 2 
uncertainty. 


4. Similar to the gradient in Sutton et al. ( |1995 1 and |Wright 


et al. (1996), CO and SO have stronger emission and higher 


abundances in the plateau of both the high- and low-velocity 
outflows, which has been suggested to be due to shock 
enhancement in these regions. Although molecules in gen¬ 
eral have around 1-2 magnitudes lower abundances in our 
low-velocity outflow than the values reported in the NE-SW 
platea u from previous interferometric study ( | Wright et al.| 
1996| ), their abundances in our high-velocity outflow are 


comparable to the single dish results in the NW-SE plateau 
( |Sutton et al.|1995| ). 


The estimated column densities of COMs (CH 3 CH 2 CN, 
CH 3 COCH 3 , CH 3 OCH 3 , and HCOOCH 3 ) in HC are f ound 
to be lower than the values in Eriedel & Snyder ( 2008| ) by a 
factor of 10 , because their observations have a smaller syn¬ 
thetic beam. However, HCOOCH 3 has comparable column 
densities to |Eavre et ( |2011a| ) in mm2 (ME4-5), mm3 a 
(ME3), and mm3b (M El), as does CH 3 OCH 3 compared to 
IBrouillet et al.| ( |2013| ) at all the central condensations (HC, 
mm2, mm3a, and mm3b), and the abundance of CH 3 COCH 3 
compared to |Peng et~n^ ( 2Q13| ) at HC and mm2 (Ace-3). 


To sum up, the abundance of most species from our study 
are comparable to previous single-dish studies especially in the 
central condensations. Therefore, molecular abundance estima¬ 
tion from single-dish observations for the farther sources (e.g., 
ten times farther than Orion-KL) should be trustworthy at least 
in the source centre after optical depth correction. 


5.2. Chemistry in the hot molecular core 

At a spatial resolution of ~ 1200 AU, our data resolves 
Orion-KL into several condensations. The specific intensities 
(Jy beam“^) of the dust continuum peaks exhibit a gradient 
from the central condensations to the extended structures, and 
the gas temperatures exhibit the same gradient: The HC is 
the strongest continuum substructure, because it has both the 
highest gas temperature and column density; mm2, mm3 a, and 
mm3b are slightly fainter and less dense; while NE and SR have 
weaker continuum and line emission than the former sources, 
but stronger emission than OEIN(S). 


Differentiation in the spatial distribution of S-, N-, and 
0 -bearing species clearly shows that chemistry in the central 
condensations of Orion-KL is not homogeneous. Simple N- and 
0 -bearing molecules show the strongest emission towards the 
central condensations. However, unlike most of the N-bearing 
species (including N-bearing COMs) which are not detected 
in SR, NE, or the outflow, simple organic molecules (H^^CO, 


CH 2 CO) and CO (^^CO, C^^O) have bright extended emission 
even towards the cooler and less dense regions (SR, OEIN, and 
OEIS). Compared to the above extreme cases, the S-bearing 
molecules are widely distributed over the central condensations, 
NE and (or) SR. A possible reason for this is that the S-bearing 
molecules are mainly formed via warm gas-phase chemistry or 
shock chem istry (e.g., |Pineau des Eorets et al.||199^ Chamley 


et al.|p^97) ), and then coincide with both N- and 0-bearing 


molecules in the plateau (NE) of BN/KL (e.g., jChandler & 
Wood| |1997[ jSchreyer & Henning |1999| 


see discussion in 
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Section |53] ). 


For the COMs, it has been widely suggested that “N-/0- 
chemical differentiation” can be caused by lower or higher 
temperatures during the cold pre-hot molec ular core phase and 
different paces of warm-up in these regions ([Caselli et al.|1993[ 

[Garrod et al.|20Q8| [Neill et aLpOlUfLaas et al.|2Qll| ). Besides 

confirming the HC peak of N-bearing COMs (CH 2 CHCN, 
CH 3 CH 2 CN, CH 3 CN, and CHf CN) and the mm3a (or/and 
mm3b) peak of several 0-bearing COMs (e.g., HCOOCH 3 
and CH 3 OCH 3 ), we found several other COMs which cannot 
be classi fied within these dichotomous emission peaks groups 
(also see |Guelin et~ST]|2008t [Friedel & Widicus Weaver |2012| 


[Peng et al.|2Q13| ). For example, the peak of CH 3 OH is velocity 

dependent, shifting from the HC (3.6 kms“^) through mm3a 
(6 kms“^) to mm3b (7.2 kms“^); the peak of CH 3 CH 2 OH is in 
between the HC and mm3 a. In contrast, CH 3 COCH 3 shares a 
morphology with N-bearing molecules - peaking towards the 
HC with a V-shaped distribution. This similarity in morpholo¬ 
gies may indicate a link between the HC-peaking N-bearing 
species and the mm3 a (or/and mm3b)-peaking species, therefore 
their formation pathways likely involve NH 3 (ammonia) or 
another major N-bearing molecule (e.g., [Chen & Woon||2011| ), 
or need similar physical conditions approp riate to their pro - 
duction and/or sublimation, such as shocks ( [Peng et 3T]|2013 ). 
Moreover, although CH 3 CHO is only tentatively detected, it 
peaks in a similar location to HCOOCH 3 . Furthermore, like 
CH 2 CO, it exhibits extended emission from the NE to the SE 
tail, suggesting these three species may share a common or at 
least a related formation path. 


The favoured model of hot molecular core chemistry so far is 
the gas-grain model, where frozen radicals, partly produced by 
photolysis, become mobil e during the warm-up phase, and their 
reacti ons produce COMs ( jMarkwick et al.||2000t [Garrod et al.j 
[2008[ ). Three phases are usually modelled, the cold, warm-up, 
and hot molecular core stages. During these stages, differ¬ 
ent species are defined as zeroth-generation (e.g., CH 3 OH), 
first-generation (e.g., HCOOCH 3 X and second-generation 
organic molecules (see [Herbst & van Dishoeck|2009[ ). In such a 
scenario, since the critical densities are fairly similar for N- and 
0-bearing species, the distribution diversity of species in Orion- 
KL may also be explained by different initial compositions or 
different chemical “ages” (evolutionary stage, the time that the 
gas phase species has had to evolve from an initial composition 
after being released from ice mantles, [Friedel & Snyder|2008| ). 
If that is the case, COMs like CH 3 OCH 3 /CH 3 OH can act as 
chemical clocks, implying that the different substructures in 
BN/KL could have different chemical ages (e.g.,[Charnley et al.[ 
[TW/ll^^am et al.|2004] ). 


In addition to grain-surface and gas-phase chemical process¬ 
ing, it is also possible that varying physical conditions, including 
temperature, density, and kinematics cause the differentiation in 
the spatial distributions of molecular line emission seen in the 
Orion-KL region. Probably due to the difference in the central 
protostars, the warm-up paces of different su bstructures vary. 
The modelling result from [Garrod et al. ( [2008 ) (especially their 
Eigure 8 ) is an example. In a particular substructure during the 
warm-up process (e.g., from 100 K to 200 K), the abundances 
of some molecules increase (e.g., HNCO, CH 3 CN, CH 3 OH) 
because of their evaporation from the grain surface to the gas 
phase, while the abundances of some large molecules cannot 
be synthesized on grains any longer because their precursor 


species have desorbed to the gas phase. They decrease (e.g., 
HCOOCH 3 , CH 3 OCH 3 ) or start to decrease (CH 3 OH). Another 
possibility is that these species are destroyed by ion-molecule 
chemistry (e.g., by protonation via XH^ molecules), followed 
by dissociative recombination into multiple small fragments. In 
this scenario, the different emission peaks of different species 
can be explained by the higher temperature in HC (120-160 K) 
than in mm3 a and mm3b (80-100 K), meaning that N-bearing 
species are still evaporating in HC but formation of molecules 
such as HCOOCH 3 is stopped or it is dissociating there. 


5.3. Chemistry in the outfiow 


Outflows from young stars can destroy grain material and/or 
liberate ices when they impact the surrounding cold matter, 
meaning that chemical rea ctions can be driven by the assoc iated 
shocks even in cold gas ( [Herbst & van Dishoeck[[200^ . Eor 
years, the outflows in Orion-KL have been the main subject 
of chemical studies of this region ( [Wright et al. 1996 and 
references therein). Although COMs are less abundant or not 
detected in these outflows, the chemistry of simple molecules 
especially shock tracers can provide an alternative way to 
investigate the gas structure of the outflows. 


Two outflows have been reported in Orion-KL, almost 
perpendicular to each other. In addition, along the NE-SW 
direction, a large-scale “integral-shaped” filament has been 


detected aligned with the low-velocity outflow ( 

Wright et al. 

\TWl |Lis et al.||1998l |Di Francesco et al.||2008 

). To confirm 


whether the NE clump is associated with the filament or the 
outflow, we selected several widespread species, and use a 
Gaussian function to fit their line profile extracted from NE 
(Eigure [^. Broad line wings (C^^O) and a second component 
(HNCO, SO 2 , and ^^SO) with Vpeak ~ 20-23 kms-i indicate 
that line emission in NE contains significant contribution from 
the outflow and is part of the plateau. 


To study the chemistry in the outflows, we image the blue- 
and red-shifted gas of several shoc k-tracing molecules along 
the outflow direction (Eigure [A4[ ). To avoid contamination 
from blending, we image as broad a velocity range as possible 
(labelled in each panel), though, -20 ~ -r30 kms“^ is the 
limit for most lines other than SO, CO, and ^^CO. Because of 
the range limit and the limits of the 52" primary beam, the 
high-velocity outflow (30- 100 kms“^) cannot be imaged in 
its entirety. However, the low-velocity outflow (~ 18 kms“^) 
with a NE-SW elongated structure is clearly detected, showing 
bipolar structure in the intensity-integrated map of H^^CO. 
In addition, intensity-integrated maps of such molecules as 
OCS, 34 s 02 , so, HNCO, CH 2 CO, CO, and ^^CO all have a 
“butterfly” morphology, with blue-shifted radial velocities of 
the moving gas towards the NW and the red-shifted velocities 
in the SE. This morphology is exactly the same as found in 


SiO ([Wright et al. 

1995 

, [Plambeck et al. 20091 

Zapata et al. 

[20121 [Niederhofer 

et al. 

[ 2012 [) and H 2 O maser 

observations 


morphology is due to the bipolar lobes of the high-velocity 
outflow or to the expanding low-velocity outflow is still not 
clear. Possibilities are discussed (e.g., two bipolar outflows from 
a processing binary, ballistic ejecta from a rotating disk, a singl e 
expanding wide-angle NE-SW outflow; Zapata et al.[[2Q12) . 
The asymmetry shape of the lobes (e.g., “U” shape blue-shifted 
and half “U” red-shifted lobes in the map of SO) is said to 
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come from some density inhomogeneities in the surroundings 
( Zapata et al.||2Q12|) or the interaction between outflow and the 
environment |Niederhofer et al.|2012| ). 

Therefore, it can be suggested that a picture of explosive 
NW-SE outflow explains the excitation temperature gradient 
from the central condensation s to the extended outflows in 
Orion-KL ( Zapata et al.||2Q12| ): the decelerated bullets can be 
absorbed by the BN/KL cloud , and their kinetic e nergy can be 
transferred to thermal energy ( |Peng et ^|2012b ); the closer a 
gas parcel is to the explosion centre (e.g., the closer encounter 
of BN, source I and maybe also source N, suggested by proper 
motio n studies from Rodriguez et al.|2005 Gomez et al.|2005 


2008| and 3D modelling from |Nissen k al.|2Q12| ), the higher the 


thermal energy available to affect the chemical properties of the 
gas in the HC and mm2. 

Of all the shock-related molecules, S-bearing species are of 
special interest, not only because their abundance may be en- 


hanced by shocks ([Martin-Pintado et al.|1992|[Pineau des Forets 

et al. 

19931 [Chemin et al.[[1994 

[Bachiller & Perez Gutierrez 

1997 

8 chilke et al.[[ 19971 Cham] 

ey et al.[[1997nBachiller et al. 

2001 

Viti et al.|2004l[Garrod|201 

13[), but also because they have 


been proposed as potential chemical clocks to date outflows 
(and hence their protostellar driving source; e.g., |Chamley et al.] 
1997t[^kelam et al.|2004[ ). The angular resolution of our study 


does not allow us to investigate the post/pre-shocked chemistry. 
In spite of that, we do And that some S-bearing molecules 
(e.g., SO 2 , ^"^802, and SO) have almost the same abundances 
in NE as in the central sources, which are higher than those in 
OEIN(S), even though we assume the temperatures in NE and 
OEIN(S) are equally low. The enhancement of their abundance 
may indicate stronger shocks in this direction than the rest 
of Orion-KL. On the other hand, although we cannot tell the 
precise ages of different substructures, we can compare the 
variations in abundance ratios SO 2 /SO in our observations to 
the gas-grain modelling result. To simplify the case, we assume 
that the temperatures in all substructures are 100 K and that 
all substructures have the same density ~ 10 ^cm“^; and we 
estimate the value of SO 2 (lines are possibl y optically thick) 
by us ing ^"^802 and a ratio of ^^ 8/^'^8 = 22 ( [Wflson & Rood 
1994[). Comparing the qbseryed 8 O 2 / 8 O ratio to modelling of 


Wakelam et al. (2004 
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in all substructures (EigurepT|), we And 
that the HC has the highest 8 O 2 / 8 O ratio, which suggests it is 
the most eyolyed substructure. Also, mm2, mm3a, and mm3b 
haye lower ratios and are younger. The lowest ratios are in the 
outflows, which indicate the age of the low-yelocity outflow 
should be 1 x 10 "^ yearK while the high-yelocity outflow 
should be ^ 1 X 10 ^ years 

The shocks could also be responsible for the obseryed 
chemical differentiation in COMs (e.g., by stripping ices from 
grains). The CR (together with mm3 a and mm3b) is thought 
to be the place where icy grain s are released into the gas 
phase by an outflow from the HC ( [Blake et ar]|1987[ [Liu et al.j 
|2002| ). Therefore, based on the reported association between the 
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In 


Wakelam et al. 


(2004), SO/H 2 ratio is sensitive to the grain com¬ 


ponents. However, SO 2 /SO ratios from different grain component mod¬ 
els are similar at a particular time (yr), so we compare this ratio to the 
observations. 

NE is part of the filamentary structure, so the age of the low- 
velocity outflow is implied from the upper limit ratio in SR. Due to the 
non-detection of SO 2 (^“^802) in OFIN(S), the age of the high-velocity 
outflow is implied from the upper limit ratio in OF IN. 



Log time (yr) 



Fig. 11: Observed 8 O 2 / 8 O ratio compared to th e model predic¬ 
tion. T he upper panel is adopted from Figure 5 of jWakelam et al. 
(2004) (100 K, 10^cm“^), where different line types are from 
different grain component models. The lower panel is from our 
observed 8 O 2 / 8 O ratio by assuming (1) all substruc tures have 
similar dens ity 10^cm“^ and (2) NSO 2 /N 34 SO 2 “ (Wilson & 


Rood |1994| ). The black solid line comes from the assumption 


that all substructures have a temperature close to 100 K, and 
the red dotted line comes from different temperatures derived in 
8 ection |4 .1 j (labelled in red). 


emi ssions of HCOQCH 3 and the 2.12 jim vibrationally excited 
H 2 ( Lacombe et al.|20()^ , it has been suggested that shocks be¬ 
tween the I 8 M and high-velocity gas “bullets” from the outflow 
( Zapata et al. 2009[|Bally et al.|201 1| ) explains the produc tion of 


COMs, whic h has the same morphology as HCOOCH 3 ( jFavre 
et al.||20lTb] ). In a similar scenario, N-bearing COMs may be 


released to the gas phase north-east of the HC (e.g., HC 3 N, 
CH 2 CHCN), where the low-velocity outflow is impacting the 
ambient dusty material, and the weak emission /nondetection 
of these molecules along the high-velocity outflow may be due 
to lack of dusty mate rial in the NW- 8 E direction ( jFriedel & 
jWidicus Weaver|2Q12| ). 


5.4. Search for COMs 

The chemistry of Orion-KL is rich, not only in high-density 
tracers and shock tracers, but also with numerous transitions of 
COMs, which are less abundant than the other tracers and which 
have a less well understood formation mechanism. 
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Three possibilities have been discussed for the formation of 
the COMs: (1) ion-molecule chemistry involving large radicals 
( [Ehrenfreund & Charnley|[^00| ); (2) formation on dust grains 
and then released to the g as at higher temperatures in the later 
stages of star formation ( [Garrod et al.||20Q8} |Herbst & van| 
Dishoeck||2009| ); (3) a hybrid approach that COMs are ejected 
into the gas via efficient reactive desorption, after a sequence 
of gas-phase reactions between precursor species formed on 


grain surfaces (|Vasyunin & Herbst 

2013|). In our dataset, 

we search for several COMs (Table 

2 ) by using the synthetic 


fitting programme described in Section [3. 2.2 1 HCyN, C 6 H, and 
CH 3 CHO are tentatively detected because of missing predicted 
lines. For the remaining species, if not taking the blending 
problem into account, all the predicted transitions are robustly 
detected. In addition, we find that saturated molecules are likely 
to be more abundant than unsaturated molecules (e.g., HC 3 N, 
C 6 H, HC 7 N). This indicates that these saturated species may be 
formed on dust grain surfaces by subsequent hydrogenation of 
their precursor species (such as CO, O, N or their more complex 
combinations). 


Unfortunately, we lack the sensitivity to detect HCOOH 
(fo rmic acid), which has been previously reported in Orion-KL 
by Liu et al^ ( |2002| ). Line blending also prevents us from ro¬ 
bustly detecting CH 2 OHCHO (glycolaldehyde) and (CH 20 H )2 
(ethylene glycol) in our data. We also fail to det ect other prebi - 
otic molecules, e.g., NH 2 CH 2 COOH (glycine; |Garrod|[2013| ), 
CH 3 COOH (acetic acid), NH 2 CH 2 CN (amin o acetonitrile; 
Belloche et al.|| 2 QQ^ , CH 3 CONH 2 (acetamide; [van Dishoeck] 


2009| ), and the ring molecules. More sensitive observations 


(e.g., with ALMA) are required to determine whether these 
molecules exist in the Orion-KL environment. 


CH 3 CHO and long carbon chains like C 6 H and HC 7 N. 

3. The Orion KL region exhibits clear chemical differentiation 
on the scales probed by the SMA-30m observations. The 
emission from S-bearing molecules, CO isotopologues, 
and simple organic molecules are extended, covering not 
only the central condensations (HC, mm2, mm3 a, and 
mm3b), but also the NE clump, the SR, and even the outflow 
regions. In comparison, the distributions of COMs are more 
concentrated and hard to group. The segregated emission 
peaks at either the HC (e.g., N-bearing) or CR (together with 
mm3a and mm3b, e.g., HCOOCH 3 , CH 3 OCH 3 ) indicates 
different formation pathways for 0-bearing and N-bearing 
molecules. However, some COMs peak in-between the 
HC and mm3 a (e.g., CH 3 CH 2 OH), or have dual velocity 
dependent peaks (CH 3 OH), or share the same morphology 
as N-bearing species (CH 3 COCH 3 ), which may indicate 
the linking of some formation pathways between 0 -bearing 
and N-bearing molecules. Moreover, line widths of the 
N-bearing, S-bearing species, HCOOCH 3 and CH 3 OH vary, 
and so do the derived rotation temperatures, indicating they 
trace different gas. 

4. By studying the chemistry within outflows of Orion-KL, we 
found higher abundances of such shock tracers as ^'^S 02 , 
SO 2 , and SO in the plateau of the low-velocity outflow 
(NE), indicating stronger shocks in this direction than the 
rest of Orion-KL. 

5. Molecular rotation temperatures and abundances show large 
gradience between central condensations and the outflow 
regions, indicating difference between hot molecular core 
and shock chemistry. 


6. Conclusions 

Based on the analysis of combined SMA and 30 m observations 
of Orion-KL, we have drawn the following conclusions: 


1. From interferometric observations with the SMA, Orion-KL 
has been resolved into several continuum substructures 
with different peak flux densities, at the angular resolution 
of 3"(1 200 AU at a distance of 414pc). In addition to the 
previous well-known compact continuum substructures like 
the HC and mm2, our data resolves the clump south-west 
of HC into two cores: mm3 a and mm3b. In addition, we 
confirmed the detection of strong continuum emissions 
(> 5 (t) from two resolved clumps (SR and NE), which lie 
along the axis of both the low-velocity outflow and the 
large-scale dusty filament. 

2. For the first time, we have combined interferometric and 
single-dish observations of Orion-KL, to correct missing 
short spacings, and achieved high sensitivity to all spatial 
scales. Covering a field of view of 52", we mapped the 
distribution of all the identified species. In particular, we 
obtained the first maps of molecules such ^^S 02 , O^^CS, 
HNCO, H^^CO, i^CO, and CH 2 CO with both high reso¬ 
lution and sensitivity to spatial information on all scales. 
We also identified several transitions of low abundance 
(< 10“^ with respect to H 2 ) COMs, such as CH 3 COCH 3 
and CH 3 CH 2 OH, as well as made the tentative detection of 
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Table Al. Identified emission lines from the SMA-30m combined dataset. 


Freq. 

(GHz) 

Mol. 

EulkB 

(K) 

Note 

Ereq. 

(GHz) 

Mol. (candidates) 

EulkB 

(K) 

Note 

218.903 

OCS(v=o)(18^ 17) 

100 






218.966 

HCOOCH3(v=1)(18i2,6 ^ 17i2,5)E 

384 






218.981 

HNCO(,=o)(10i,io ^ 9 i, 9) 

101 






219.068 

HCOOCH3(v=i)(18i7,2 ^ 17i7,i)E 

481 






219.079 

HCOOCH3(v=1)(283,25 ^ 282,26)E 

434 






219.090 

HCOOCH3(v=0)(347,28 ^ 345,29)E 

387 

* 

219.090 

HCOOCH3(y=o)(347 28 ^ 345 29)A 

387 


219.154 

HCOOCH3(v=i)(18ii,7 ^ 17ii,6)E 

369 

* 

219.153 

HCOOCH3(y=i)(104,6 ^ 93,6)E 

230 


219.174 

HC3N(,,=i)(24^23) 

452 

t 





219.195 

HCOOCH3(v=1)(18i6,3 ^ 17i6,2)E 

459 

* 

219.196 

HCOOCH3(y=o) (366,30 ^ 366,3 i)E 

429 


219.220 

CH3COCH3(v=0)(21i,20 ^ 202,19)AE 

122 

* 

219.220 

CH3COCH3(y=o)(2l2,20 ^ 20i,i9)AE 

122 






219.220 

CH3COCH3(y=0)(2l2,20 ^ 202,19)EA 

122 






219.220 

CH3COCH3(y=o)(21i,20 ^ 20i,i9)EA 

122 


219.242 

CH3COCH3(v=0)(21i,20 ^ 202,19)EE 

122 

* 

219.242 

CH3COCH3(y=o)(2l2,20 ^ 20i,i9)EE 

122 






219.242 

CH3COCH3(y=0)(2l2,20 ^ 202,19)EE 

122 






219.220 

CH3COCH3(y=o)(21i,20 ^ 20i,i9)EE 

122 


219.264 

CH3COCH3(v=0)(21i,20 ^ 202,19)AA 

122 

* 

219.264 

CH3COCH3(y=o)(2l2,20 ^ 20i,i9)AA 

122 


219.276 

SO2(v=0)(227J5 ^ 236,18) 

352 






219.311 

CH3COCH3(v=0)(129,4 ^ ll8,3)AA 

66 






219.325 

UL 


7 





219.331 

HCOOCH3(v=1)(18i5,4 ^ 17i5,3))E 

438 

* 

219.328 

HCOOCH3(y=i)(369,28 ^ 368,29)E 

635 


219.355 

^"^SO2(v=0)(lll,ll ^ 10o,1o) 

60 

t 





219.401 

CH2CHCN(v=0)(233,20 ^ 223,19) 

145 






219.412 

HCOOCH3(v=i)(18io,8 ^ 17io,7)E 

355 






219.417 

HCOOCH3(y=0) (305,26 ^ 304,27)E 

291 






219.463 

CH3CH2CN(222,21 ^21i,20) 

112 






219.479 

HCOOCH3(y=i)(18i4,5 ^ 17i4,4)E 

419 






219.484 

HCOOCH3(y=0)(305,26 ^ 304,27)A 

291 






219.506 

CH3CH2CN(y=0)(242,22 ^ 232,2l) 

136 






219.547 

HNCO(y=o)(104,6 ^ 94,5) 

750 

* 

219.547 

HNCO(y=o)(104,7 ^ 94,6) 

750 


219.560 

Ci^O(2 ^ 1) 

16 

t 





219.566 

HCOOCH3(y=i)(18i5,4 ^ 17i5,3)A 

438 

* 

219.566 

HCOOCH3(y=i)(18i5,3 ^ 17i5,2)A 

438 


219.568 

HCOOCH3(y=i)(18i4,5 ^ 17i4,4)A 

419 

* 

219.568 

HCOOCH3(y=i)(18i4,4 ^ 17i4,3)A 

419 


219.571 

HCOOCH3(y=i)(18i6,3 ^ 17i6,2)A 

458 

* 

219.571 

HCOOCH3(y=i)(18i6,2 ^ 17i6,i)A 

458 


219.579 

HCOOCH 3 (y=o)( 289,19 ^ 288 , 2 o)A 

294 

* 

219.579 

HCOOCH3(y=i)(18i7,2 ^ 17i7,i)A 

481 






219.579 

HCOOCH3(y=i)(18i7,i ^ 17 i7,o)A 

481 






219.578 

HCOOCH3(y=i)(299,2i ^ 298,2i)E 

498 


219.584 

HCOOCH3(y=i)(18i3,6 ^ 17i3,5)A 

401 

* 

219.584 

HCOOCH3(y=i)(18i3,5 ^ 17i3,4)A 

401 






219.586 

HCOOCH3(y=o(309,22 ^ 308,23)A 

330 


219.592 

HCOOCH 3 (y=o( 289 ,i 9 ^ 288,2o)E 

294 






219.600 

HCOOCH3(y=o(309,22 ^ 308,23 )E 

330 






219.607 

HCOOCH3(y=o(305,26 ^ 303,27)E 

291 

* 

219.610 

CH3COCH3(y=0)(3 3 8,25 ^ 338,26)EE 

368 






219.610 

CH3COCH3(y=0)(3 39,25 ^ 337,26)EE 

368 


219.622 

HCOOCH3(y=i)(18i2,7 ^ 17i2,6)A 

384 

* 

219.622 

HCOOCH3(y=i)(18i2,6 ^ 17 i2,5)A 

384 


219.642 

HCOOCH3(y=i)(18i3,6 ^ 17i3,5)E 

401 






219.657 

HNCO(103 8 ^937) 

447 

* 

219.657 

HNCO(103,7 ^ 93,6) 

447 


219.696 

HCOOCH3(y=i)(18ii,8 ^ 17 ii,7)A 

369 

* 

219.696 

HCOOCH3(y=i)(18u,7 ^ 17 u,6)A 

369 


219.705 

HCOOCH3(y=i)(184,i5 ^ 174,14)A 

299 






219.734 

HNCO(y=o)(102,9 ^ 92,8) 

231 






219.737 

HNCO(y=0)(102,8 ^ 92,7) 

231 






219.764 

HCOOCH3(y=i)(189,9 ^ 179,8)E 

342 






219.798 

HNCO(y=o) (160,10 ^ 60,9) 

58 

t 





219.822 

HCOOCH3(y=i)(18io,9 ^ 17io,8)A 

355 

* 

219.822 

HCOOCH3(y=i)(18io,8 ^ 17io,7)A 

355 


219.827 

HCOOCH3(y=i)(18i2,7 ^ 17i2,6)E 

384 






219.909 

hPC0(3i,2^2i,i) 

33 

t 





219.949 

SO(y= 0)(65 ^ 54 ) 

35 

t 





219.984 

CH3OH(yj=0) (253,22 ^ 244,2o)E 

802 






219.994 

CH3OH(yj=0)(23549 ^ 226 ,i7)E 

776 






220.030 

HCOOCH3(y=i)(189,io ^ 179,9)A 

342 

* 

220.030 

HCOOCH3(y=i)(189,9 ^ 179,8)A 

342 


220.043 

HCOOCH3(y=i)(18ii,8 ^ 17 ii,7)E 

368 






220.078 

CH30H(yj=o)(8o,8 ^ 2i,6)E 

97 






220.154 

t-CH3CH20H(243 22 ^ 242 23) 

263 






220.167 

HCOOCH3(y=0)(174,i3 ^ 164,12)E 

103 






220.178 

CH2CO(11i,ii^10i,io) 

76 

t 





220.190 

HCOOCH3(y=o)(174,i3 ^ 164 , 12 )A 

103 






220.258 

HCOOCH3(y=i)(l88,10 ^ 178,9)E 

330 






220.296 

CH13CN(129 ^ II9) 

646 






220.307 

HCOOCH3(y=i)(18io,9 ^ 17io,8)E 

354 






220.324 

CH3CN(y=o)(12io^llio) 

782 






220.355 

CH3COCH3(y=o)(22o,22 ^ 21o,2i) 

124 

AE,EA* 

220.355 

CH3COCH3(y=o)(22i,22^21i,2i) 

124 

AE,EA 





220.355 

C6H(J = 159/2-157/2,^=3/2, l=e) 

425 

7 

220.362 

CH3COCH3(y=o)(22o,22 ^ 21 i, 2 i)EE 

124 

* 

220.362 

CH3COCH3(y=o)(22i,22 ^ 21i,2i)EE 

124 






220.362 

CH3COCH3(y=o)(22i,22 ^ 21o,2i)EE 

124 






220.362 

CH3COCH3(y=o)(22o,22 ^ 21o,2i)EE 

124 
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220.368 

HCOOCH3(v=i)(188,ii ^ 178,io)A 

331 

* 

220.370 

HCOOCH3(y=i)(l88,10 ^ 178,9)A 

331 





220.368 

CH3COCH3(y=o)(22i,22 ^ 21 i,2i)AA 

124 





220.368 

CH3COCH3(y=o)(22o,22 ^ 21o,2i)AA 

124 





220.368 

CHi3CN(128^1l8) 

526 

220.399 

13C0(2 ^ 1) 

16 

t 




220.404 

CH3 CN(,=o )(129 ^ II9) 

646 





220.416 

HCOOCH3(v=1)(183,16 ^ 172,15)A 

293 





220.431 

CH13CN(127 ^ II7) 

417 

* 

220.433 

HCOOCH3(y=0)(355,28 ^ 335,29)A 

357 

220.446 

CH3CHO(133,io^ 132,11V = 0)E 

105 





220.465 

CH3COCH3(y=o)(llii,i ^ 10io,o)AE 

63 

* 

220.466 

CH3COCH3(y=o)(llii,o ^ 10 io,i)AE 

63 

220.476 

CH3CN(y=o)(128 ^ llg) 

526 





220.486 

CH13CN(126 ^ ll6) 

326 

* 

220.487 

C6H(J = 159/2-157/2, n=3/2„ \=f) 

426 

220.525 

HCOOCH3(v=i)(104,6 ^ 93,7)A 

231 





220.533 

CHf CN(125 ^ II5) 

247 





220.539 

CH3CN(y=o)(127 ^ II7) 

419 





220.561 

CH2CHCN(24i,24 ^ 23i,23) 

134 





220.570 

CH13CN(124 ^ II4) 

183 





220.585 

HNCO(y=o)(10i,9 ^ 9 i,8) 

lOl 





220.594 

CH3CN(y=o)(126 ^ ll6) 

326 





220.621 

CHj3CN(122^1l2) 

97 





220.641 

CH3CN(y=o)(125 ^ II5) 

248 





220.661 

CH3CH2CN(y=o)(252,24 ^ 242,23) 

143 





220.679 

CH3CN(y=o)(124 ^ II4) 

183 





220.709 

CH3CN(y=o)(123 ^ II3) 

133 





220.730 

CH3CN(y=o)(122 ^ II2) 

98 





220.743 

CH3CN(y=o)(12i^lli) 

76 





220.747 

CH3CN(y=o)(12o ^ llo) 

69 





220.764 

CH3COCH3(v=o)(11ii,o ^ 10io,o)EE 

63 





220.786 

HCOOCH3(y=0)(283,25 ^ 283,26)E 

248 





220.812 

HCOOCH3(y=o)(183,i6 ^ 172,15)E 

105 





220.815 

HCOOCH3(v=0)(183,16 ^ 172,15)A 

105 





220.820 

HCOOCH3(v=0)(242,23 ^ 24i,24)A 

169 





220.848 

CH30CH3(244,20 ^ 235,19) 

297 

A(E)A(E) 




220.866 

HCOOCH3(v=0)(283,25 ^ 283,26)A 

248 

* 

220.866 

HCOOCH3(y=i)(192,i7 ^ 183 ,i6)A 

303 

220.889 

HCOOCH3(v=0)(18i7,i ^ 17i7,o)A 

293 

* 

220.889 

HCOOCH3(y=o)(18i7,2 ^ 17 i7,i)A 

293 

220.893 

CH30CH3(234,20 ^ 233,21) 

274 

A(E)A(E) 




220.901 

HCOOCH3(y=o)(18i7,i ^ 17 i7,o)E 

293 





220.913 

HCOOCH3(v=i)(187,i2 ^ 177,ii)A 

320 

* 

220.910 

HCOOCH3(y=o)(18i7,2 ^ 17 i7,i)E 

293 

220.926 

HCOOCH3(v=o)(18i6,2 ^ 17i6,i)A 

271 

* 

220.926 

HCOOCH3(y=o)(18i6,3 ^ 17 i6,2)A 

271 

220.935 

HCOOCH3(y=o)(18i6,2 ^ 17 i6,i)E 

271 





220.946 

HCOOCH3(v=i)(187,ii ^ 177,io)A 

321 

* 

220.947 

HCOOCH3(y=o)(18i6,3 ^ 17 i6,2)E 

271 

220.967 

HC7N(y=o)(J= 196-195) 

1045 

7 




220.978 

HCOOCH3(v=0)(18i5,4 ^ 17i5,3)A 

250 

* 

220.978 

HCOOCH3(y=o)(18i5,3 ^ 17 i5,2)A 

250 

220.985 

HCOOCH3(y=i)(187,ii ^ 177 ,io)E 

320 

* 

220.984 

HCOOCH3(y=o)(18i5,3 ^ 17 i5,2)E 

250 

220.998 

HCOOCH3(y=o)(18i5,4 ^ 17 i5,3)E 

250 

* 

220.999 

g-CH3CH20H(13o,i3 ^ 12o,i2,vt = 0-1) 

136 

221.048 

HCOOCH3(v=0)(18i4,5 ^ 17i4,4)A 

231 

* 

221.048 

HCOOCH3(y=o)(18i4,4 ^ 17 i4,3)A 

231 





221.048 

HCOOCH3(y=o)(18i4,4 ^ 17 i4,3)E 

231 

221.067 

HCOOCH3(v=0)(18i4,5 ^ 17i4,4)A 

231 





221.076 

HCOOCH3(v=0)(299,21 ^298,22)A 

312 





221.086 

HCOOCH3(y=o)(299,2i ^ 298,22)E 

312 





221.111 

HCOOCH3(y=i)(188,ii ^ 178 ,io)E 

330 





221.115 

^"^SO2(v=0)(222,20 ^ 22 i,2i) 

248 





221.124 

CH2CHCN(y=o)(23i,22 ^ 22i,2i) 

130 





221.140 

HCOOCH3(y=o)(18i3,5 ^ 17 i3,4)E 

213 

* 

221.141 

HCOOCH3(y=o)(18i3,5 ^ 17 i3,4)A 

213 





221.141 

HCOOCH3(y=o)(18i3,6 ^ 17 i3,5)A 

213 

229.056 

CH3COCH3(v=o)(22i,2i ^2l2,2o)EE 

133 

* 

229.056 

CH3COCH3(y=o)(22i,2i ^21i,2o)EE 

133 





229.056 

CH3COCH3(y=o)(222,2i ^ 2l2,2o)EE 

133 





229.056 

CH3COCH3(y=o)(222,2i ^21i,2o)EE 

133 





229.058 

CH3COCH3(y=o)(149,6 ^ 138,5)EE 

85 

229.078 

CH3COCH3(v=o)(22i,2i ^ 21 2,20)AA 

133 

* 

229.078 

CH3COCH3(y=o)(222,2i ^21i,2o)AA 

133 

229.087 

CH2CHCN(y=o)(243,2i ^233,20) 

156 





229.117 

HCOOCH3(y=i)(193,i7 ^ 182,16)E 

303 





229.127 

CH3COCH3(v=0)(12io, 2 ^ ll9,2)EE 

68 





229.134 

CH3CH0(93,7 ^ 92,7) 

62 





229.203 

UL 


7 




229.224 

HCOOCH3(v=o)(239,i4 ^ 238,i5)A 

217 





229.259 

HCOOCH3(y=o)(239,i4 ^ 238,i5)E 

217 

* 

229.258 

CH3CHO(103,7 ^ 102,8) 

71 

229.265 

CH3CH2CN(y=o)(262,25 ^ 255,24) 

154 





229.320 

HCOOCH3(y=o)(239,i5 ^ 238,i6)E 

217 





229.348 

SO2(v=0) (115,7 ^ 124,8) 

122 

t 




229.389 

HCOOCH3(y=o)(239,i5 ^238,i6)A 

217 





229.405 

HCOOCH3(y=o)(183,i5 ^ 173,14)E 

111 





229.420 

HCOOCH3(y=o)(183,i5 ^ 173,14)A 

111 





229.474 

HCOOCH3(y=o)(203,i7 ^ 194 ,i6)E 

134 





229.491 

t-CH3CH20H(175,i2 ^ 174,13) 

160 

* 

229.498 

g-CH3CH20H(4111,31 ^ 40i2,29,vt = 1-0) 

925 





229.498 

g-CH3CH20H(4111,30 ^ 40i2,28,vt = 1-0) 

925 

229.505 

HCOOCH3(y=o)(203,i7 ^ 194 ,i6)A 

134 
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229.539 

HCOOCH3(v=0)(3l5,27 ^ 3l3,28)E 

309 






229.589 

CH3OH(vt=0)(154,ii ^ 163J3)E 

374 






229.607 

HCOOCH3(v=0)(3l5,27 ^ 3l3,28)A 

309 






229.648 

CH2CHCN(v=0)(25i,25 ^24i,24) 

146 






229.661 

HCOOCH3(v=i)(259,i6 ^ 258,i7)A 

432 






229.759 

CH30H(,^=0)(8-i,8^7o,7)E 

89 






229.858 

^SO2(v=0)(42,2 ^ 3i,3) 

19 






229.864 

CH3OH(yj=0)(195,15 ^ 204J6)A + + 

578 






229.888 

^'^SO2(v=0)(173,15 ^ 18o,18) 

162 






229.939 

CH3OH(yj=0)(195,14 ^ 204J7)A-- 

578 






229.987 

^''SO2(v=0)(426,36 ^ 417,35) 

921 






230.027 

CH3OH(y^=0)(3-2,2^4_i,4)E 

40 






230.109 

C 6 H(J = 165/2-163/2, ^=1/2, \=f) 

482 






230.140 

CH30CH3(254,22 ^ 253,22)EA 

319 

A(E)A(E)$ 





230.170 

CH3COCH3(y=o)(23o,23 ^ 22 i, 22 )AE 

135 

* 

230.170 

CH3COCH3(y=o)(23i,23 ^ 220 , 22 )AE 

135 






230.170 

CH3COCH3(y=o)(23o,23 ^ 22 o, 22 )EA 

135 






230.170 

CH3COCH3(y=o)(23i,23 ^ 22 i, 22 )EA 

135 


230.177 

CH3COCH3(y=o)(23o,23 ^ 22 i, 22 )EE 

135 

* 

230.177 

CH3COCH3(y=o)(23i,23 ^ 22o,22)EE 

135 






230.177 

CH3COCH3(y=o)(23o,23 ^ 22 o, 22 )EE 

135 






230.177 

CH3COCH3(y=o)(23i,23 ^ 22 i, 22 )EE 

135 


230.183 

CH3COCH3(y=o)(23o,23 ^ 22 i, 22 )AA 

135 


230.183 

CH3COCH3(y=o)(23i,23 ^ 220 , 22 )AA 

135 


230.234 

CH30CH3(172,15 ^ 163,14) 

148 

A(E)A(E) 

230.234 

CoHCJ = 165/2-163/2, n=l/2, l=e) 

483 

7 

230.283 

CH3COCH3(y=0)(168,8 ^ 159,7)EE 

107 






230.293 

CH3OH(yj=0) (222,20 ^ 21_3J9)E 

609 






230.318 

0‘3CS(19^ 18) 

111 

t 





230.369 

CH3OH(y,=0)(224,i8 ^ 2 l 5 ,i 7 )E 

683 






230.377 

HCOOCH 3 (y=o)( 229 j 4 ^ 228 ,i 5 )A 

203 






230.436 

HCOOCH3(y=o)(25i,24 ^ 25i,25)E 

182 

* 

230.436 

HCOOCH3(y=o)(25i,24 ^ 25o,25)E 

182 


230.441 

HCOOCH3(y=o)(252,24 ^ 25i,25)E 

182 

* 

230.441 

HCOOCH3(y=o) (252,24 ^ 25o,25)E 

182 


230.481 

CH3OCH3(108,2^1l7,5)EE 

140 

* 

230.474 

CH30CH3(108,3 ^ ll7,4)EA 

140 






230.476 

CH30CH3(108,3 ^ ll7,4)EE 

140 






230.478 

CH30CH3(108 3 ^ ll7 4)AA 

140 






230.478 

CH30CH3(108;2 ^ ll7!5)AA 

140 






230.479 

CH30CH3(108 3 ^ ll 7 4)AE 

140 






230.479 

CH 3 OCH 3 ( 108 ; 2 ^ 1 l 7 ! 5 )AE 

140 






230.483 

CH3OCH3(108,2^1l7,5)EA 

140 


230.488 

CH2CHCN(y=o)(24i,23 ^23i,22) 

141 






230.505 

CH30CH3(264,23 ^ 255,20) 

343 

A(E)A(E) 





230.538 

CO(y=0)(2 ^ 1) 

17 

t 





230.673 

g-CH3CH20H(132,ii ^ 122 , 10 ,vt = 0 - 0 ) 

139 






230.697 

g-CH 3 CH 20 H (83 5 ^ 827 ,n = 1 - 0 ) 

103 






230.739 

CH2CHCN(yLo)(25o,25 ^ 24o,24) 

146 






230.794 

g-CH3CH20H(65,i ^ 54 , 1 ,Vt = 0-1) 

105 

* 

230.794 

g-CH 3 CH 20 H( 65,2 ^ 54 , 2 ,n = 0-1) 

105 


230.844 

HCOOCH3(y=i)(19i7,2 ^ 18i7,i)E 

493 






230.852 

HCOOCH3(y=i)(19i6,3 ^ 18i6,2)E 

470 






230.879 

HCOOCH3(y=i)(184,i4 ^ 174,13)A 

301 






230.889 

HCOOCH3(y=i)(19i5,4 ^ 18i5,3)E 

450 






230.933 

^^SO 2 (v= 0 )( 54,2 ^ 63 , 3 ) 

52 

* 

230.936 

HCOOCH3(y=0)(293,26 ^ 292 , 27 )E 

265 


230.954 

t-CH 3 CH 20 H( 165 ,ii ^ 164 , 12 ) 

145 






230.960 

HCOOCH3(y=i)(19i4,5 ^ 18i4,4)E 

430 






230.991 

t-CH3CH20H(14o,i4 ^ 13 i , i 3 ) 

86 


230.992 

g-CH 3 CH 20 H( 65 ,i ^ 64 , 2 ,n = 1-1) 

110 






230.993 

g-CH 3 CH 20 H( 65,2 ^ 64 , 3 , Vt = 1-1) 

110 


231.019 

HCOOCH3(y=0)(124,9 ^ ll 3 , 8 )E 

57 

* 

231.020 

HCOOCH3(y=0)(293,26 ^ 292 , 27 )A 

264 


231.046 

HCOOCH3(y=0)(124,9 ^ ll 3 , 8 )A 

57 






231.061 

OCS(y=o)(19^18) 

111 

t 





231.102 

HC 7 N(J = 205-204) 

1142 






231.188 

HCOOCH 3 (y=o)( 2 l 9 ,i 3 ^ 2 l 8 ,i 4 )E 

190 






231.199 

HCOOCH 3 (y=o)( 2 l 9 ,i 2 ^ 2 l 8 ,i 3 )A 

190 

* 

231.199 

HCOOCH3(y=0)(2l9,i2 ^ 2 l 8 ,i 3 )E 

190 


231.221 

i3CS(y=o)(5^4) 

33 

t 





231.231 

HCOOCH3(y=i)(19i2,7 ^ 18i2,6)E 

396 


231.231 

g-CH3CHO(128,5 ^ ll8,4E) 

216 



Note. 1. Due to the frequency resolution (0.8125 MHz, ~ 1.2kms“^), lines with broad FWHM width can be attributed to different species, so lines with 

stronger CDMS/JPL intensity are listed here in the left column with and the potential blended weaker transitions are listed in the right column. 

2. Tentative detections and unidentified lines are marked with 

3. Lines with “f” are imaged in Figure|^ with “ij:” are imaged in Figur ej^ 

4. “A(E)A(E)” represents 4 types of transitions are possible: AA, EE, AETEA. 
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Table A2: Laboratory parameters of CH 3 CN and CH^^CN and the intensity integrated over the width of each line f TB(v)dv, which is used for rotation diagrams in Figure 
Uncertainties on the measured intensities are typically < 10% (written as the subscript), as determined from Gaussian fitting. 
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Table A3: The laboratory parameters and the intensity integrated over the width of each line f TB(r')dL', which is used for rotation diagrams in FigureUncertainties on the 
measured intensities are typically < 10% (written as the subscript), as determined from Gaussian or HFS fitting. 
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Table A3: (continued) 
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Integrated intensity f Tb (v)dv (K km s ^) 
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Table A6: Column densities and abundances for N-bearing molecules from different substructures in Orion-KL denoted in Figure [^and[^ 
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Table A7: Column densities and abundances for S-bearing molecules from different substructures in Orion-KL denoted in Figure]^ 
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Specific Intensity (Jy/beam) 
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229 229.5 230 




Rest Frequency (GHz) 


Fig. Al: Spectrum from the HC, with identified strong lines marked. From 220.5 GHz to 220.8 GHz, CH 3 CN and CH^^CN are 
labelle d in red and blue, respectively. COMs containing only C, H, O elements and other tentatively detections are labelled in 
Figure [A^ 
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CH 3 OH @ mm3a (Trot=160 K; Nmol = 2.5 x 10*'' cm-^; Vis, = 7.6 kms"'; AV = 5 kms"') 



Trot = 102 K; Nmol = 1.0E+17 cm^2; AV = 3 km/s 






219.05 219.1 

Frequency (GHz) 


219.15 219.2 


219.4 219.45 

Frequency (GHz) 




220.03 220.04 

Frequency (GHz) 


220.17 220.18 220.19 

Frequency (GHz) 


Fig. A2: Sy nthe tic spectrum fitting of COMs towards their emission peaks (Visr of each species are corrected according to Vpeak 
from Figure [A3] ). The black histograms are the observed spectra, overlaid with a best fit model spectrum for given uniform FWHM 
line width (aV ) (fit from Figure [Xs] ), Trot (assumed from an adjacent substructure), and Nmol (adjusted based on Nt in Tables [A5^ , 
|A6^ , and j A7P listed on top of each panel. Blue line mark the central frequency of each transition. Uncertainties of the fitting to 
observations come from the assumption that under LTE, the candidate molecule lines are optically thin, and each species has the 
same line width at different transitions. Referenced molecular data are from CDMS/JPL ( Sanchez-MongepOTT Palau et al.|2Qll I. 
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220.25 220.3 220.35 

Frequency (GHz) 



220.75 220.8 220.85 220.9 220.95 

Frequency (GHz) 



Frequency (GHz) 




221 221.05 221.1 221.15 

Frequency (GHz) 



230.84 230.86 230.88 230.9 230.92 230.94 

Frequency (GHz) 



230.96 230.98 231 231.02 231.04 

Frequency (GHz) 

CH 3 COCH 3 @ SMAl (Trot=126 K; N^oi = 1-5 x 10*® cm-2; V^r = 5.8 kms 



231.18 231.2 231.22 

Frequency (GHz) 


-1; AV = 3kms“i) 



Fig.A2: (continued) 
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230.15 230.2 230.25 230.3 

Frequency (GHz) 

CH 3 OCH 3 @ mm3b (Trot=92 K; N^oi = 2.0x 10^^ cm-^; Visr = 5.6 kms-^; 



AV = 4 kms 




CH 3 CH 2 OH @ near SMAl & mm3a (Trot=70 K; N^oi = 1-5 x lO^^ cm-^; Visr = 8.4 kms-^; AV = 2.5 kms-^) 







Fig.A2: (continued) 
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CH 2 CHCN @ North to hotcore (Trot=155 K; N^oi = 5.0x10'^ cm-^; = 4.9 kms-i; AV = 7 kras"') 









Frequency (GHz) 


CH 3 CH 0 @ mmSb (Trot=102 K; N^oi = 1-3 x 10*® cm-^; Vis, = 8.0 kms-*; AV = 5 kms"*) 




CgH @ SMAl (Trot=126 K; N^oi = 1.0x 10*® cm-^; = 4.0 kms"*; AV = 4 kms"*) 




HC 7 N @ North to hotcore (Trot=155 K; N^oi = l.Ox 10*® cm-^; Vis^ = 6.3 kms"*; AV = 5 kms"*) 




220.95 220.96 220.97 220.98 

Frequency (GHz) 


231.09 231.1 231.11 

Frequency (GHz) 


Fig.A2: (continued) 
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(Jy/beam) 
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Fig. A3: Line profiles of identified species from each continuum peak. In measuring the main Vpeak and AV of each line at each 
substructure, lines in green and notes in blue (in the form of “value(error)”) are given by the Gaussian fittings using Gildas software. 
Multi-Gaussian fittings are used to CH 3 CN, CH^^CN, HCOOCH 3 (y=o), and CH 3 COCH 3 due to blending, and to SO 2 , ^"^802, and 
SO due to the second velocity component (the second velocity component of HNCO measured from multi-Gaussian fittings are 
given in Figure |7]). Single Gaussian fittings are used to the rest lines, while lines whose Vpeak cannot be measured by Gaussian 
fittings owing to strong self-absorption are marked with and non-detections are not marked. Two red vertical lines in each 
panel mark the velocity range over which the emission from each species is integrated to make the distribution maps shown in 
Figures]^ and 1^ 
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(Jy/beam) 



Fig. A3: (continued) 
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Fig. A3: (continued) 
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CH3OCH3 CH3COCH3 CH3CH2OH CH3CHO CgH 

(254,22-253,22) ( 23 o,23-221,22) ( 14 o. 14 - 15 i,, 3 ) (1 53,10 -1 52 ,i 1) (J= 1 65 / 2-1 63 / 2 ) 
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Fig. A3: (continued) 


HNCO HC 3 N(^ 7 ^i) CH 3 CN CH 3 CN CH 2 CHCN CH 3 CH 2 CN HC 7 N 

(10o,,o-9o,9) (24-23) (122-1 1 2 ) (122-11 2 ) (25o,25-24o,24) ( 262 . 25 - 252 , 24 ) (J = 205-204) 
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Fig. A3: (continued) 
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20 0 -20 20 0 -20 20 0 -20 

Aa (") 

Fig. A4: Line wing flux integration maps from the combined SMA-30 m data at 1.3 mm. The grey maps show the intensity integrated 
ranging through the main emission of their linewidth in Figure |A3[ Red and blue contours show the red-shifted and blue-shifted 
gas, with the intensity integrated as labelled. Yellow lines sketch the outflow directions. Black crosses denote the hotcore, SMAl, 
Source I, and Source N; white crosses denote the BN object and CR. 
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Feng et al.: Resolving the chemical substructure of Orion-KL 
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Fig.AS: Channel maps of the second velocity component (Vpeak ~ 22 km s of the strongest unblended ^"^802 and SO2 line, as 
well as of 80(65 ^ 54). 
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Feng et al.: Resolving the chemical substructure of Orion-KL 
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CH 3 OH ( 195 ,i 4 - 204 ,i 7 )A— 229.939 GHz 
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Fig.A6: Channel maps of unblended CH 3 OH lines show multiple velocity-dependent emission peaks: 2.4-3.6 kms“^ towards 
SMAl, ~ 4.8-6 kms“^ towards HC, 7.2-9.6 kms“^ towards HC/mm3a/mm3b, 10.8 kms“^ towards mm3a/mm3b/mm2 and 
> 12 kms“^ towards mm2. Green lines sketch the outflow directions. Black crosses denote the hotcore, SMAl, Source I, and Source 
N; white crosses denote the BN object and the CR. 


51 















